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ABSTRACT 

We present results of X-ray observations of a sample of 15 clusters selected via their imprint on 
the cosmic microwave background (CMB) from the thermal Sunyaev-Zel'dovich (SZ) effect. These 
clusters are a subset of the first SZ-selected cluster catalog, obtained from observations of 178 deg^ 
of sky surveyed by the South Pole Telescope. Using X-ray observations with Chandra and XMM- 
Newton, we estimate the temperature, Tx, and mass, Mg, of the intracluster medium (ICM) within 
r.soo for each cluster. From these, we calculate Yx = Mg Tx and estimate the total cluster mass using 
a M500 — Yx scaling relation measured from previous X-ray studies. The integrated Comptonization, 
Yszi is derived from the SZ measurements, using additional information from the X-ray measured 
gas density profiles and a universal temperature profile. We calculate scaling relations between the 
X-ray and SZ observables, and find results generally consistent with other measurements and the 
expectations from simple self-similar behavior. Specifically, we fit a Ysz — Yx relation and find a 
normalization of 0.82 ± 0.07, marginally consistent with the predicted ratio of Ysz/Yx = 0.91 ± 0.01 
that would be expected from the density and temperature models used in this work. Using the 
Yx derived mass estimates, we fit a Ysz — M^oo relation and find a slope consistent with the self- 
similar expectation of Ysz oc M^^^ with a normalization consistent with predictions from other X-ray 
studies. We compare the X-ray mass estimates to previously published SZ mass estimates derived 
from cosmological simulations of the SPT survey. We find that the SZ mass estimates are lower by 
a factor of 0.89 ± 0.06, which is within the ^^15% systematic uncertainty quoted for the simulation- 
based SZ masses. Overall, the X-ray measurements confirm that the SZ-selected sample consists of 
very massive systems which exhibit general properties consistent with X-ray selected samples, even 
allowing for the broad redshift range (0.29 < z < 1.08) of the sample. 
Subject headings: galaxies: clusters — X-rays: galaxies: clusters 
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Large area cluster surveys extending to high redshift 
can be used to study the evolution of the abundance 
of galaxy clusters, thereby delivering precise constraints 
on the amount and nature of the dark energy (Wang 
& Steinhardt 1998; Haiman et al. 2001). The accuracy 
with which the observed mass proxy can be linked to 
the true cluster mass fundamentally limits the cosmo- 
logical constraints from the survey. In particular, a red- 
shift dependent bias on a survey's cluster mass estimates 
could mimic a time-evolving dark energy, so this system- 
atic must be understood and constrained. Clusters are 
known to evolve — through mergers, galaxy and star 
formation, and variable contributions from active galac- 
tic nuclei (AGN) — and these effects will influence the 
evolution of any cluster observable at some level. 

Most of the baryonic mass in clusters is in the form 
of an intra-cluster gas that can be heated to several keV 
as it virializes in a massive cluster's gravitational po- 
tential well. This gas is visible in the X-ray band via 
thermal bremsstrahlung and also from its distortion of 
the Cosmic Microwave Background (CMB) from inverse 
Compton scattering, otherwise known as the Sunyaev- 
Zel'dovich (SZ) effect (Sunyaev & Zeldovich 1972). The 
lowest scatter cluster observables that scale with cluster 
mass, M, are likely those most closely related to the gas 
pressure (e.g., Kravtsov et al. 2006), and hence related to 
the total energy of the gas. The SZ intensity is propor- 
tional to the Comptonization, the line-of-sight integral 
of the gas pressure. Hence, the SZ signal integrated over 
the cluster's extent, Ygz, measures the total pressure in 
the cluster. An X-ray analogue, Yx, can be constructed 
from the product of the cluster's total gas mass and the 
X-ray spectroscopic temperature (Kravtsov et al. 2006). 

The intrinsic scatter in the M — Yx relation has been 
found to be smaller than the measurement error in X-ray 
studies conducted to date (Vikhlinin et al. 2006, 2009a; 
Sun et al. 2009), and has been used to place interest- 
ing constraints on the dark energy equation of state with 
only a relatively small sample of clusters (Vikhlinin et 
al. 2009b). The slope and scatter of the relationship be- 
tween Ysz and the cluster mass, M, has been studied ex- 
tensively in simulations and is expected to be relatively 
insensitive to non-gravitational physics (Nagai 2006), the 
dynamical state of clusters (Jeltema et al. 2008), and the 
presence of cool cores (Motl et al. 2005). The expected 
close correlation between cluster mass and Ysz, as well as 
the redshift independence of the SZ brightness, strongly 
motivates using SZ cluster surveys for cosmological stud- 
ies (Carlstrom et al. 2002). 

Recently, there has been significant progress in measur- 
ing the SZ-signal from clusters. High resolution imaging 
has been obtained for single objects (e.g., Nord et al. 
2009; Mason et al. 2010) and ICM profiles have been mea- 
sured for moderately sized samples (e.g., Mroczkowski et 
al. 2009; Plagge et al. 2010). The first clusters discovered 
with a blind SZ survey were reported in Staniszewski et 
al. (2009), and showed the capability of the SZ-signal as 
a cluster finder. From larger samples, scaling relations 
have been measured between the SZ-signal and mass es- 
timates from both X-ray and weak lensing measurements 
(e.g., Bonamente et al. 2008; Marrone et al. 2009; Mefin 
et al. 2010). Notably, Melin et al. (2010) measured a Ysz- 
M relation from binned WMAP fiuxes at the location of 
known X-ray selected clusters combined with X-ray lu- 



minosity based mass estimates. They found a Ysz-M re- 
lation with a normalization and slope that matched the 
X-ray prediction, however similar analyses have found 
confiicting results (Komatsu et al. 2010). More detailed 
observations comparing SZ and X-ray measurements will 
improve our understanding of the gas pressure in clusters. 

The first SZ-selected catalog of clusters was presented 
in Vanderlinde et al. (2010) (hereafter VIO), and in- 
cluded the first meaningful cosmological constraints from 
an SZ cluster survey. In VIO, the sample of 21 clus- 
ters had a median redshift of z = 0.74 and was pre- 
dicted to be 100% complete above a mass threshold of 
M500 « 3 X lO^^h-'^MQ at z > 0.6. The cosmologi- 
cal constraints from VIO were limited by uncertainties in 
the cluster mass calibration. This calibration relied on 
the dark matter simulations of Shaw et al. (2009), with 
gas physics based on the models in Bode et al. (2007), to 
link the SZ-significance to cluster mass. This introduced 
a ^ 15% systematic uncertainty in the mass calibration 
due to uncertainties in the pressure normalization of the 
simulations. Therefore, an important first step to im- 
prove the cosmological constraints of VIO is to tie the 
SZ observables to observationally calibrated X-ray scal- 
ing relations, such as those in Vikhlinin et al. (2009a). 
The VIO cluster sample is also unique because of its high 
median redshift and SZ-selection. X-ray observations of 
these clusters could allow additional constraints on the 
redshift evolution of X-ray scaling relations, which typi- 
cally have been studied from X-ray selected samples con- 
centrated at z < 0.6 (e.g., Maughan 2007; Vikhlinin et 
al. 2009a; Mantz et al. 2009). 

In this work, we present results from X-ray observa- 
tions of a subset of 15 clusters with the highest SZ- 
significance from VIO. We report on the X-ray observ- 
ables of this sample, and use X-ray scaling relations from 
Vikhlinin et al. (2009a) to estimate each cluster's mass. 
For each cluster, we also note details of the cluster's X- 
ray morphology and the identification of X-ray sources 
with objects in other catalogs. We use the X-ray mea- 
sured gas density profiles to improve the SZ estimates of 
integrated Comptonization, Ygz- Finally we construct 
X-ray and SZ scaling relations, specifically the Ysz-Yx 
and Ysz-Mc,QQ relations, and compare these relations to 
expectations and other results. 

The paper is organized as follows: §2 describes the clus- 
ter sample and the SZ, optical and X-ray observations. §3 
discusses the X-ray data analysis, the estimation of gas 
mass and temperature while §4 describes the estimate 
of gravitational mass. In §5 we discuss the analysis of 
the SPT data and deprojection of the SZ measurements. 
In §6 the X-ray and SZ scaling relations are investigated 
and we conclude with a discussion in §7. The properties 
of individual clusters are discussed in the Appendix. 

In all calculations we have assumed a 
WMAP7-^BAO+iJo ACDM cosmology (Komatsu 
et al. 2010) with VIm = 0.272, I^a = 0.728 and 
Hq = 70.2 km s~^ Mpc"^ with distance measurements 
from Baryon Acoustic Oscillations (BAO) in the distri- 
bution of galaxies (Percival et al. 2010) and the Hubble 
constant {Hq) measurement from Riess et al. (2009). 
Everywhere, we define Afsoo as the mass inside r5oo, 
within which the matter density is 500 times the critical 
density at the cluster redshift, Pcrit{z) — 3H'^{z)/8ttG, 
where H{z) is the Hubble parameter. 
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2. CLUSTER SAMPLE AND OBSERVATIONS 

The clusters used in this work were selected from a 
SZ detection significance-limited catalog from the South 
Pole Telescope (SPT) cluster survey described in VIO. 
For all clusters, we performed follow-up optical imag- 
ing to identify galaxy cluster counterparts and measure 
redshifts (High et al. 2010). Sixteen cluster candidates 
with SZ detection significance above 5.4 were selected for 
X-ray follow-up with Chandra and XMM-Newton. One 
candidate was later discovered to be a spurious detec- 
tion, and is discussed in more detail in §A. This section 
briefly describes the SZ, optical, and X-ray data sets. 

2.1. SPT SZ Observations and the Cluster Sample 

The 10-m diameter South Pole Telescope is a millime- 
ter wavelength telescope located at the South Pole. Its 
primary science goal is to conduct a ~2000 deg^ survey 
to find clusters of galaxies via measurements of the SZ 
effect. The receiver consists of a 960 element bolometer 
array that is sensitive in three frequency bands centered 
at 95, 150, and 220 GHz. In this work, we only use ob- 
servations at 150 GHz, the SPT frequency band with the 
most SZ sensitivity. Details of the telescope and receiver 
can be found in Padin et al. (2008) and Carlstrom et al. 
(2009). 

In 2008, the SPT surveyed two '-.^100 deg^ regions. 
These two fields are approximately square on the sky, 
and centered at right ascension (R.A.) 5^30™, declina- 
tion (decl.) -55° and R.A. 23'*30", decl. -55°. The 
data from the first of these fields was used to report the 
first SZ-discovered clusters (Staniszewski et al. 2009), to 
measure source counts of extragalactic mm-wavelength 
emitting objects (Vieira et al. 2009), and to measure 
small scale temperature anisotropics due to the SZ effect 
from unresolved clusters and emission from point sources 
(Lueker et al. 2009; Hall et al. 2009). In VIO, the data 
from both fields was used to report the first SZ detection 
significance-limited catalog of 22 cluster candidates. 

The SZ observations, data processing, and mapmaking 
used in this work are described in detail in VIO, and only 
a brief overview is provided here. Each field was scanned 
in azimuth at a constant velocity, with the scans stepped 
in elevation. With the velocity and elevation step used, 
it takes ^2 hours to cover an entire 100 deg^ field. This 
scan is then repeated several hundred times to decrease 
the noise in the co-added map. Each detector timestream 
is filtered to remove both long timescale drifts and sky 
signal that is spatially correlated across the focal plane. 
The filtering effectively acts as a high pass filter; the filter 
used for the SZ cluster analysis removes signal on spatial 
scales larger than ^ 0.5°. The data is combined to make 
a map by reconstructing the pointing for each detector 
and then averaging the data from all the detectors using 
inverse-variance weighting. The 150 GHz maps were cal- 
ibrated to an accuracy of 3.6% by direct comparison to 
the WMAP 5-year maps (Lueker et al. 2009). The final 
map of each field has a sensitivity limit of 18 /iK-arcmin. 

Cluster candidates were identified in the SPT maps 
by using a matched spatial filter technique (Haehnelt & 
Tegmark 1996; Mclin et al. 2006). Here we summarize 
the method and results as used in VIO. The SPT maps 
consist of several sources of signal, including primary 
CMB anisotropy, unresolved point sources, and SZ signal 



from clusters. To identify cluster candidates, the SPT 
maps are filtered in Fourier space to give more weight- 
ing to signals matching the expected spatial scales of 
clusters. Twelve different spatial filters were constructed 
using spherical /3-modcls with (3 fixed to 1 and core radii 
evenly spaced between 0.25' and 3.0'. With these, the 
SPT maps were filtered and cluster candidates were iden- 
tified as decrements in the filtered map. The significance 
of a candidate was quantified by their signal relative to 
the standard deviation in the filtered map, or signal-to- 
noise. For a given cluster candidate, the highest signal- 
to-noise across all filter scales was defined as ^. To avoid 
spurious identifications from ringing around spatially fil- 
tered bright point sources, a 4' radius region around all 
5(7 positive sources was masked before spatial filtering. 
The total sky area used in both fields after masking was 
178 deg^. Further details of the SZ analysis relevant to 
the results presented here are given in §5. 

The sixteen cluster candidates with the largest ^ were 
selected for X-ray follow-up with the Chandra and XMM- 
Ncwton X-ray telescopes. The original candidate list was 
based on a similar, but earlier version of the list that ap- 
peared in VIO. This change re-ordered the list somewhat, 
such that sixteen of the seventeen most significant can- 
didates from VIO had X-ray follow-up. In Table 1, we 
give the position and ^ from VIO for the fifteen confirmed 
clusters with X-ray observations. One of the candidates 
with X-ray follow-up, SPT-CL J2343-5521, is not fisted 
because it is very likely a false detection. Its X-ray ob- 
servation is discussed further in the Appendix. 

2.2. Optical Imaging and Spectroscopy 

Optical counterparts for the clusters selected by SZ de- 
tection significance were identified and photometric red- 
shifts were measured via a combination of imaging from 
the Blanco Cosmology Survey (BCS) and targeted ob- 
servations using the Magellan telescopes. For a subset of 
clusters, spectroscopic redshifts were also obtained. Fur- 
ther details of the optical data and analysis can be found 
in High et al. (2010), and are briefly described below. 

The BCS is an optical survey of two ^50 deg^ flclds 
that he inside the two 2008 SPT fields described in §2.1. 
BCS used the Mosaic-II wide field imager on the Blanco 
4-m telescope at the Cerro Tololo Inter- American Obser- 
vatory in Chile, in Chile. The BCS obtained contiguous 
deep optical imaging in the griz bands across their sur- 
vey fields, and these data have been processed using a 
development version of the DES data management sys- 
tem (Ngeow et al. 2006; Mohr et al. 2008) and then used 
in the study of the galaxy population and in the redshift 
estimation of the first SPT survey fields (Ngeow et al. 
2009; Staniszewski et al. 2009; High et al. 2010; Brodwin 
et al. 2010; Zenteno et al. 2010). These data arc pub- 
licly available through the NOAO Survey Program and 
have been used by other groups to study clusters in the 
SPT survey region (Mcnanteau et al. 2009; Mcnanteau 
& Hughes 2009; Mcnanteau et al. 2010; Mclnnes et al. 
2009; Suhada et al. 2010). For the nine clusters out- 
side the BCS fields, optical imaging in the griz bands 
was obtained with the Inamori Magellan Areal Camera 
and Spectrograph (IMACS; Dressier et al. 2006) on the 
Magellan Baade 6.5-m telescope in Las Campanas, Chile. 
Five cluster candidates in the BCS fields were also re- 
observed with Magellan. In contrast to BCS, the Mag- 
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ellan observations were performed adaptively, where the 
candidates were observed in '-^100 sec increments until 
the galaxy cluster was detected. 

Optical counterparts of each SPT candidate were iden- 
tified by searching for red sequence objects within a 2' 
radius of the SPT candidate location. A cluster was iden- 
tified through an excess of red sequence objects relative 
to the background, and the photometric redshift was es- 
timated by fitting a red sequence model. The redshift 
uncertainty varies over the sample, however they are typ- 
ically Az ^ 0.03 and can be as large as Az = 0.10 for 
clusters at z ~ 1 (see High et al. 2010, for details). 

For 8 of 16 cluster candidates, spectroscopic measure- 
ments were obtained using the Low Dispersion Survey 
Spectrograph (LDSS3) on the Magellan Clay 6.5-m tele- 
scope. These observations are described in High et al. 
(2010). For one candidate, SPT-CL J0546-5345, multi- 
slit spectroscopy was done with IMACS on the Magel- 
lan Baade 6.5-m telescope (Brodwin et al. 2010). The 
spectroscopic targets were chosen to span the redshift 
range of the sample to help calibrate the photometric 
redshifts. One cluster, SPT-CL J0516-5430, has a pre- 
vious spectroscopic redshift from the REFLEX cluster 
survey (Bohringcr et al. 2004). In Table 1, we give the 
photometric redshifts for each of the 15 clusters included 
in this work, and the spectroscopic redshift where avail- 
able. 

2.3. Chandra & XMM-Newton observations 

As described in §2.1, the sixteen cluster candidates 
with the highest detection significance, ^, were selected 
for an X-ray follow-up program. This program was split 
between several observing cycles and proposals between 
the Chandra and XMM-Newton X-ray satellites. 

In the original planning of X-ray observations, the ob- 
serving time required for each cluster was estimated from 
preliminary predictions of the SZ significance-to-mass re- 
lation and a mass-luminosity scaling. The uncertainty 
in the relations does not allow an accurate estimate of 
the exposure needed for a required number of photons. 
The estimates were informed by ROSAT fluxes where 
available. For clusters scheduled for multiple observa- 
tions with Chandra, the integration time was modified 
accordingly once a flux had been measured. 

To date, the Chandra program consists of two GTO 
programs in AO-9 (295 ks, total), two GTO programs in 
AO-11 (340 ks, total) and a GO program in AO-11 (310 
ks). When completed, the Chandra program will have 
collected at least 1500 cluster photons within 0.5r5oo and 
in the 0.5-7.0 keV energy band for each of 12 clusters in 
the sample. This limit was chosen to enable measuring of 
the ICM temperature, Tx, to 15% accuracy. For one very 
faint object, SPT-CL J0553-5005, the data wiU not be 
sufficient to measure the temperature to this accuracy in 
the current program. The Chandra ACIS-I count-rate in 
the 0.5-7.0 kcV band for this object is only 0.005 ct s~^. 

Four of the candidates have been observed with XMM- 
Newton, one of which is an archival observation (SPT- 
CL J0516-5430, also RXCJ0516. 7-5430); the other three 
are from a 2008 program (SPT-CL J0559-5249, SPT-CL 
J2332-5358 and SPT-CL J2337-5942). AH of these ob- 
servations have more than the required 1500 counts. For 
two of the clusters, there exists both Chandra and XMM- 
Newton data (SPT-CL J0516-5430 and SPT-CL J2337- 



5942) and these are all analyzed in this work. We chose 
to re-observe these two clusters with Chandra to better 
identify and remove X-ray point-sources in the analysis. 
The low redshift and high mass of the clusters made these 
observations possible with a comparably small amount of 
observing time. 

While we find consistent results from the analysis of 
data from the two different satellites, we choose to use the 
Chandra data for the two clusters with data from both 
instruments since both the high spatial resolution and 
stable background of Chandra are desirable. For the two 
clusters where only XMM-Newton data were available, 
we derive our results from these data using the methods 
as described below. 

X-ray measurements of SPT-CL J2332-5358 and SPT- 
CL J2342-5411 were also reported by Suhada et al. (2010) 
from the XMM-BCS survey. These measurements were 
not as deep as the observations discussed in our work, 
however we find X-ray observables consistent with their 
measurements. 

Currently, there are eight out of the fifteen clusters for 
which the required 1500 source counts have not yet been 
collected with either instrument. However, we include 
these in the analysis since the existing data provide useful 
constraints. 

3. X-RAY DATA ANALYSIS 
3.1. Data reduction 

The exposure times and resulting source counts from 
the Chandra and XMM-Newton observations are listed 
in Table 1. We also list the cluster coordinates and SZ 
signal-to-noise as described in VIO as well as the opti- 
cal redshifts and the galactic absorbing hydrogen column 
from the Leiden- Argentine-Bonn survey (Kalberla et al. 
2005). Source counts are quoted within 0.5r5oo which 
is estimated from the Yx parameter and the M^qq-Yx 
relation (see §3.5). 

The Chandra data were reduced using CIAO 4.1 and 
CALDB 4.1.3. AU data were taken with the ACIS-I nom- 
inal aimpoint in VFAINT telemetry mode and additional 
screening to reject particle background was applied. To 
remove periods of fiaring background we extracted point- 
source subtracted lightcurves in the 0.3-12 keV band and 
filtered these using a 3a threshold. XMM-iVcwton data 
were reduced using SAS 9.0 and reprocessed. Source- 
free light curves were generated in hard (MOS:10-12 keV, 
pn: 12-14 keV) and soft (MOS:0.3-10 keV, pn:0.3-10 keV) 
energy bands separately and a 3(t cut was applied to re- 
move periods of high background. 

3.2. Data analysis methods 

For the X-ray observables, Tx, Mg and Yx we use scal- 
ing relations from Vikhlinin et al. (2009a) iteratively to 
determine the rsoo radius where the observable is mea- 
sured. 

The rsoo radii used for the measurement of Tx were 
estimated using the M-Tx relation in Vikhlinin et al. 
(2009a) from samples of local clusters with deep expo- 
sures for which hydrostatic masses could be determined: 

M5oo = (3.02 ±0.11) 10^^ h-^ Mq 

/ ,rp \ (1.53±0.08) 
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TABLE 1 

Chandra and XMM-Newton observations 





R.A'' 


Dec.'' 


£^ 
s 


photo-^ 


spec-"" 


exposure 


Source 




rl V a 

— 


Merger? 
















counts 


10^" cm^^ 


t 




SPT-CL JOOOO-5748 


0.25 


-57.807 


5.48 


0.74 


- 


28.1 


1451 


1.37 


0.16 




SPT-CL J0509-5342 


77.336 


-53.705 


6.61 


0.47 


0.4626 


27.3 


2441 


1.46 


0.13 


/ 


SPT-CL J0516-5430 


79.148 


-54.506 


9.42 


0.25 


0.2952S 


8.4 


2136 


2.05 


0.15 


/ 


SPT-CL J0528-5300 


82.017 


-53.0 


5.45 


0.75 


0.7648 


36.5 


356 


3.23 


0.20 




SPT-CL J0533-5005 


83.398 


-50.092 


5.59 


0.83 


0.8810 


41.5 


201 


2.95 


0.53 


/ 


SPT-CL J0546-5345 


86.654 


-53.761 


7.69 


1.16 


1.0665^ 


55.6 


1304 


6.78 


0.13 


/ 


SPT-CL J0551-5709 


87.902 


-57.156 


6.13 


0.41 


0.4230 


19.8 


876 


6.27 


0.53 


/ 


SPT-CL JO559-5249X 


89.925 


-52.826 


9.28 


0.66 


0.6112 


42.0 


2006 


5.08 


0.42 


/ 


SPT-CL J2331-5051 


352.958 


-50.864 


8.04 


0.55 


0.5707 


34.2 


2428 


1.12 


0.23 


/ 


SPT-CL J2332-5358^'' 


353.109 


-53.976 


13.05 


0.32 




18.9 


4826 


1.28 


0.36 




SPT-CL J2337-5942 


354.354 


-59.705 


14.94 


0.77 


0.7814 


19.8 


1488 


1.45 


0.12 


/ 


SPT-CL J2341-5119 


355.299 


-51.333 


9.65 


1.03 


0.9983 


79.0 


2090 


1.21 


0.26 




SPT-CL J2342-5411 


355.69 


-54.189 


6.18 


1.08 




133.7 


1193 


1.49 


0.23 




SPT-CL J2355-5056 


358.955 


-50.937 


5.89 


0.35 




20.6 


1798 


1.27 


0.61 




SPT-CL J2359-5009 


359.921 


-50.16 


6.35 


0.76 




57.9 


713 


1.33 


0.78 


/ 



^X.MM-Newton observation. The listed exposure times and source counts are for MOSl, MOS2 and PN combined. 
''R.A. and Dec. determined from the SPT detection. 

Signal- to- noise measured in 150 GHz SPT maps as described in VIO. 
■^Listed photometric and spectroscopic redshifts are based on analysis in High et al. (2010). 
"^Based on Chandra/XMM-Newton count-rate in the 0.5-7.0 keV band, within 0.5 rsoo 
•^Hydrogen column density from the Leiden-Argentine-Bonn survey (Kalberla et al. 2005). 
^Distance between SPT detection and X-ray centroid. 

s Spectroscopic redshift from the REFLEX cluster survey (Bohringer et al. 2004). 
'^Spectroscopic redshift from (Brodwin et al. 2010). 

'The position and ^ of this cluster are modified from VIO after correcting for the effects of a point source coincident with the 
cluster, as described in 55.4. 



where 



= V(l + 2)^(1 + ^Mz) - z{2 + z)nA. 
The radius is then defined as 



'"500 = 



/ 3M: 



500 



yAiTbOOpcntiz) 



1/3 



(2) 



(3) 



Similarly, in the estimation of the gas mass, Mg, de- 
scribed below, we use the gas-fraction relation 



/3,5oo = (0.0764 ± 0.004) h~^-' 

(0.0225 ± 0.004) /i^i-MogAfis (4) 

where /g,5oo is the gas mass fraction within rsoo, /g.500 = 
Mg/M^QO, and M15 is the total mass, A/500, m units of 
10^^ Mq'^^. With this relation we exphcitly take 
into account the observed trend of /g.500 with cluster 
mass (e.g., Vikhlinin et al. 2009a). The gas mass estima- 
tion is much more dependent on the aperture radius than 
the estimate of temperature and Mg must be estimated 
iteratively to obtain a self-consistent result. Given a gas 
mass, we use equation 4 to determine the total mass and, 
hence, r5oo. 

For the estimation of the Yx = Mg x Tx parameter 
we analogously estimate the aperture radius through the 
determination of 



A/500 = (5.77 ± 0.20) 10" h^/'^MQ 

\ (0.57±0.03) 

3 X 10" A/okeV ) ^^^^ 



-2/5 



, (5) 



^® This equation contains a typo in Vikhlinin et al. (2009a), 
Table 3. 



also from Vikhlinin et al. (2009a), and iteratively deter- 
mine Yx and the total mass within rsoo (See §3.5 below). 

For each cluster, the center of the X-ray emission is 
determined using a centroid calculation. X-ray point 
sources were excluded in this process. The distance be- 
tween the X-ray centroid and the SPT position is listed 
in Table 1. The SZ and X-ray positions given are both 
centroid measurements. However, due to the different 
weighting of their signals on the ICM density and tem- 
perature, the centroids are not expected to agree for a 
cluster which is not azimuthally symmetric. For Chan- 
dra data, point sources were identified using the CIAO 
tool wavdetcct and removed from subsequent analysis. 
In the XMM-Newton data, point sources were detected 
with the SAS task edetect chain. Additionally, extended 
secondary maxima were identified and removed following 
the method prescribed in Vikhlinin et al. (2009a). Ex- 
tended substructures were included in the estimation of 
the total cluster luminosity. 

3.3. Spectral analysis 

Spectra and response files were generated from the 
Chandra data for the spectral extraction regions using 
specextract. We use two independent methods of back- 
ground subtraction in this work and compare the results. 
Due to the limited spatial extent of most of the SPT can- 
didates it is possible to use in-field background subtrac- 
tion. Background was extracted from regions near the 
source but outside of the r5oo radius where the source 
flux is a small fraction of the background flux. Since the 
nominal pointing was chosen to be offset from the clus- 
ter position it is possible to extract the background at 
a detector coordinate where the mirror effective area is 
similar to that of the source extraction region and where 
the particle-induced background is similar. This limits 
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complications related to under-subtracted X-ray back- 
ground which can be a problem when using in- field back- 
ground. As a second background subtraction method we 
also use the latest blank field backgrounds collected by 
M. Markevitch^^, re-project these on the sky to match 
our cluster datasets and normalize the exposure times to 
match the count rates in the 9.5-12 keV band. The back- 
ground spectra can then be subtracted from the same 
detector region as the source. We compare the temper- 
atures and flux estimates from these two techniques and 
find good agreement. Since the blank- field background 
subtraction provides better signal-to-noise, we use this 
method throughout as this allows for temperature esti- 
mates within the r^oo aperture for most objects. 

For the XMM-Ncwton observations we use the spectral 
analysis as described in Andersson & Madejski (2004) 
and Andersson et al. (2009). For this work we limit 
the analysis to use in-field background since we are only 
interested in obtaining one spectrum for each cluster. 
Spectra for both Chandra and XMM-Ncwton datasets 
were extracted inside r5oo, excluding the central emis- 
sion within 0.15r5oo to avoid the effects of cool cores, 
known to cause additional scatter in X-ray scaling rela- 
tions. We model the data using a MEKAL model for 
the thermal plasma with emission lines and a WABS ab- 
sorption model. For our sample, the data are not deep 
enough to measure metal abundances reliably and we fix 
the abundance to 0.3 solar. The hydrogen equivalent ab- 
sorbing column, uh, was fixed at the weighted average 
from the Leiden- Argentine-Bonn survey (Kalberla et al. 
2005). These are listed in Table 1. The best fit temper- 
atures are shown in Table 2. 

3.4. X-ray imaging 

The X-ray surface brightness is extracted in n concen- 
tric annuli defined by 



r, = I IX 



l/c 



1 ■ ■ - n, 



(6) 



where we have used r^ax = 1-5 ^500, c = 1.2 and n = 20. 
The limiting radius Vmax is chosen to be large enough 
so that confusion between background and cluster flux is 
negligible and so that the integrated gas density is not 
overestimated due to projection effects. The values of c 
and n are chosen to balance high central resolution with 
achieving a similar number of X-ray counts per radial 
bin. 

For every annular region the average exposure is cal- 
culated, taking into account bad pixels and chip gaps. 
Similarly, the average effective mirror area is calculated 
for each region taking into account the mirror vignetting. 
The radial model for the X-ray counts includes a compo- 
nent for the cluster X-ray surface brightness, a spatially 
flat unvignettcd component for the particle induced X- 
ray background and soft X-ray background plus an ad- 
ditional spatially flat vignetted component taking into 
account under-subtraction of the soft X-ray background 
flux due to mirror vignetting. This last component is 
always small. 

The cluster surface brightness is taken to be pro- 
portional to the integrated emission measure, EM = 

http:/ /cxc. harvard. cdu/contrib/maxim/acisbg 



J UeUpdV, which is a good approximation in the 0.5 — 2 
keV band for gas temperatures present in our sample. 
The inferred density distribution has been shown previ- 
ously to have little dependence on the temperature when 
determined from the surface brightness in this energy 
band (Mohr et al. 1999; Forman et al. 2007). The ra- 
dial gas density distribution is modeled using a modified 
/3-model, 



(r/re)-" 



'(l-|-r2/r-?)3'3-"/2 {l + r-r/riyh'' 



(7) 



(Vikhlinin ct al. 2006). The model accounts for the cuspy 
centers of clusters as well as the steepening of the profile 
seen at larger radii. We fix 7 = 3 in the above expression 
(see Vikhlinin et al. 2006). 

For the XMM-Newton data, the projected model is 
convolved with a model of the XMM point-spread func- 
tion (PSF), as described in Ghizzardi et al. (2001), prior 
to fitting. The model is normalized by integrating equa- 
tion 7 over a cylindrical volume with a radius of r^oo, 
truncated at 3r5oo along the line-of-sight. This is com- 
pared to the spectral normalization of a MEKAL model 
with parameters derived in §3.3 within r^oo- This way 
the model normalization, no, can be calculated using the 
angular distance Da determined from the redshift listed 
in Table 1. The gas mass Mg.500 is calculated by setting 
He = Zup and pg = mpUeA/Z where Z ~ 1.199 and 
A = 1.397 are the average nuclear charge and mass re- 
spectively for a plasma with 0.3 solar abundances (imply- 
ing a mean molecular weight per free electron = 1.165) 
assuming abundances by Anders & Grevesse (1989). The 
gas density pg is integrated over a sphere of radius r^oo 
determined iteratively for Mg and Yx estimates using 
-^500 — Mg and AI^qq — Yx relations, respectively. 

The parameter space of the model is explored via 
Markov chain Monte-Carlo iteration and the uncertain- 
ties on the parameters are estimated using the Markov 
chain posterior. 

3.4.1. Cluster images 

X-ray images in the 0.5 - 2.0 kcV band are generated 
using the filtered event files and binned into 4" x 4" pix- 
els. These are shown in Figures 5-19 alongside optical 
images of the clusters from High et al. (2010). Images 
are smoothed with a 8" Gaussian filter. The images are 
shown with SZ S/N contours from VIO and the white 
crosses show the location of the brightest cluster galaxy 
(BCG). In all cases, the SPT detection is within 1' of 
the X-ray centroid and in a majority of cases (9 of 15) it 
is within 20". Additionally, Sydney University Molonglo 
Sky Survey (SUMSS, Mauch et al. 2003) radio sources 
are marked in the images with small circles (15" radius). 
Extended X-ray structures are marked with larger cir- 
cles (30" radius) or arrows and are discussed for each 
individual cluster in the Appendix. 

3.5. ICM temperature, gas mass, luminosity and Yx 

The average spectral temperature of the intra-cluster 
medium is measured within the [0.15-l]r5oo aperture as 
described in §3.3. A new value of t^qq is then estimated 
using this temperature with equation 1 and a new spec- 
trum is extracted using the new radius. The process is 
repeated until the value of Tx has converged. For 5 of 
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the clusters, the number of source photons within [0.15- 
Ijrsoo is less than 1000 and the signal-to-noise is less than 
20. This leads to large uncertainties in the spectral fits 
and could potentially cause systematic biases in the tem- 
perature. For these objects, we instead use the estimate 
of Tx within [0.15-0.5]r5oo, where the signal-to-noisc is 
higher and extrapolate using the fitted relation between 
the two temperatures from Vikhlinin et al. (2009a): 

TV,(0.15-l)r5oo/^-X,(0.15-0.5)r5oo = 

0.9075 + 0.0625 T^,(o.i5-o.5)r™o- («) 
We also add a 10% uncertainty to these temperature es- 
timates to account for the uncertainty in this relation. 
The clusters for which this procedure was performed are 
marked in Table 2. For SPT-CL J0551-5709, a cluster 
which is spatially coincident with the local cluster Abcll 
S0552, we also perform this procedure, but with a dif- 
ferent motivation. The surface brightness to the south 
and north-west of this cluster suggests that there is a sig- 
nificant amount of fiux from cooler gas associated with 
the Abell cluster. This is also seen in the temperature 
estimates. Tx drops from 4.4 ± 0.7 keV to 3.3 ± 0.4 
keV when comparing the inner and outer apertures. We 
include SPT-CL J0551-5709 in the scaling relation fits 
using the corrected temperature and mark is with a red 
square in the scaling relation plots. 

Similarly, the gas mass within r5oo, -A^g.sooi is esti- 
mated as described in §3.4 using an initial estimate of 
7-500 from equation 1 and the spectrally derived temper- 
ature, Tx- The estimate of r5oo is then revised using 
the implied gas mass fraction and equation 4, and the 
process is repeated until the gas mass converges. 

In this work, we adopt redshift evolution of /g,500i in 
the same way as it was applied in the work of Vikhlinin 
et al. (2009a). The redshift dependence of the gas- mass 
fraction is not well constrained observationally to high- 2; 
and independent measurements of the mass are needed 
to study this further. The fg^^m-M relation is used here 
only to determine the radius rsoo within which to esti- 
mate Mg and has no impact on the Yx -based mass esti- 
mates. Wc do not quote any total masses based on the 
fg.500-M relation in this work. 

SPT-CL J0516-5430 has an unusually extended mor- 
phology compared to other objects in the sample and the 
assumption of spherical symmetry is very approximate. 
The line-of-sight dimension for this cluster and the in- 
ferred gas mass are likely to be overestimated for this 
reason (see, e.g. Nagai et al. 2007). To keep the analysis 
analogous to Vikhlinin et al. (2009a) we do not attempt 
to correct for this here. SPT-CL J0516-5430 is marked 
with a blue triangle in the scaling relation plots. 

When estimating Yx, we use equation 5 to determine 
r5oo and re-calculate both Tx and Mg in a similar way. 
Cluster simulations show that Yx exhibits less scatter 
with mass compared to Mg and Tx individually due to 
the typical anti-correlation of deviations from the mean 
of these two observables (Kravtsov et al. 2006). This 
makes the Yx mass estimator less dependent on the dy- 
namical state of clusters, as far as it can be estimated 
from simulations. The X-ray observables are listed in 
Table 2. 

For some of the clusters in the sample, the X-ray data 
are not deep enough to detect the X-ray emission with 
high significance out to v^qq. For five of the clusters, the 



signal-to-noise in our radial bins (equation 6) drops below 
3 at a radius of about 0.6 r.soo- These are the same five 
objects with low signal-to-noisc spectra described above. 
The gas density profile is then constrained primarily by 
the central surface brightness profile and extrapolated 
out to rsoo. This leads to systematic uncertainties that 
are difficult to account for as they depend on cluster mor- 
phology. In our sample we find a typical density slope 
Pg oc r°' oi a ~ —1.90 ± 0.30 at 0.6 rsoo, where we quote 
the mean and standard deviation. Varying the slope out- 
side of this radius by ±0.30 typically changes the mass 
outside 0.6 rsoo by 8%. This is also an over-estimate of 
the uncertainty since density profiles generally become 
steeper with increasing radius out to rsoo- To be con- 
servative, we add an uncertainty of 10% on the total 
gas mass for clusters with poor signal-to-noise outside 
0.6 rsoo. 

4. TOTAL MASS ESTIMATES 

Using the mass proxies Tx and Yx, we estimate the 
gravitational mass, Afsoo for the clusters using equations 
1 and 5. The mass estimates are listed in Table 3. Clus- 
ters are identified as mergers based on their morphology. 
Clusters with secondary maxima, filamentary structure 
or significant isophotal centroid shifts are classified as un- 
relaxed (e.g., Mohr et al. 1993). As discussed in Kravtsov 
ct al. (2006), the Tx based masses should be multiplied 
by a factor 1.17 for clusters identified as mergers. Fol- 
lowing Vikhlinin et al. (2009a) we correct our masses up- 
ward with this factor and add an uncertainty of 5% on 
the masses (Table 3). Our merger classification is listed 
in Table 1. 

In Tabic 3, we also give the SZ-derived mass estimates 
from VIO. In VIO, the SPT cluster survey was used to 
constrain cosmological parameters while simultaneously 
fitting a cluster detection significance-mass relation. The 
significance-mass relation had priors imposed on its nor- 
malization, slope, and redshift evolution that were moti- 
vated by simulated thermal SZ maps of the sky. These 
maps were generated from large dark matter simulations 
(Shaw et al. 2009) that used a semi-analytic gas model 
of Bode et al. (2007) which was calibrated such that the 
simulated clusters matched observed X-ray scaling rela- 
tions for low redshift {z < 0.25) clusters. The maxi- 
mum likelihood significance-mass relation was then used 
to generate mass estimates for each cluster, including 
both statistical and systematic uncertainties. These sys- 
tematic uncertainties were equivalent to a ^15% uncer- 
tainty on the mass estimate for each cluster. 

Comparing the X-ray and SZ-based mass estimates, 
we find that the SZ-derived masses are lower by a factor 
of 0.83 ± 0.06, where we use only the statistical uncer- 
tainties to quote the uncertainty on the average ratio 
from a fit to a scaling relation between the two with 
the slope fixed to 1. Two clusters are 2a inconsistent 
with unity and both have larger X-ray mass estimates; 
SPT-CL J0516-5430 and SPT-CL J0546-5345. VIO notes 
that clusters at z < 0.3 will have mass estimates biased 
low because the power-law scaling that is assumed for 
the significance-mass relation does not fully capture the 
effects of CMB-confusion on the SZ signal. This effect 
could possibly explain the relatively low SZ-inferrcd mass 
found for SPT-CL J0516-5430, which is the only cluster 
in this sample at z < 0.3. Leaving this cluster out, the 
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TABLE 2 

X-RAY OBSERVABLES Tx,500, A^s.SOO, Yx ,500 AND Lx,50o(0.5 - 2.0 keV) 



Name 



JOOOO- 


■5748'' 


0, 


.74 


J0509- 


■5342 


0, 


.4626 


J0516- 


■5430 





.2952 


J0528- 


■5300'' 


0, 


.7648 


J0533- 


■5005'' 


0, 


.8810 


J0546- 


■5345 


1, 


.0665 


J0551- 


■5709'' 


0, 


.4230 


J0559- 


■5249 


0, 


.6112 


J2331- 


■5051 


0, 


.5707 


J2332- 


■5358 


0, 


.32 


J2337- 


■5942 


0, 


.7814 


J2341- 


■5119 


0, 


.9983 


J2342- 


■5411" 


1, 


.08 


J2355- 


■5056 


0, 


.35 


J2359- 


■5009" 


0, 


.76 



^X,500 


10" a/q 


'^X,500 




-^X.500 


keV 


10" Mq keV 


lO** erg s-i 


8.6_2 3 


. 7+0.5 
^•"-0.6 


4.7± 1, 


.8 


0. i it U.4 


"0-1.1 


c £1+0.2 


4.3 ± 0, 


.8 


1 n 1 


9-8-1.2 


17.0_(,4 


15.9 ± 2, 


.4 


4.7 di 0.1 


c n + 3.5 

5-2_i.7 


Q+0.4 


1.7±0, 


.9 


i.y it U.z 


4.0_i 2 


n -,+0.5 
2-1-0.4 


0.9 ±0, 


.4 


1.2 di 0.3 


7 5+1-^ 
'-''-l.l 


7 q+0.4 
'-■^-0.3 


5.3 ± 1, 


.0 


6.4 di 0.4 


4 1+0-9 


K 1+0.3 


2.0 ±0, 


.4 


1.9 it 0.2 


'•'-0.8 


O q+0.3 


6.1 ±0, 


.8 


3.3 it 0.2 


5 9+1=* 


c 7+0.2 
^-'-0.2 


3.5 ±0, 


.6 


4.4 it 0.2 




5 fi+0-2 
0.0_Q 2 


4.4 ±0, 


.6 


3.0 it 0.1 




9 5+0-4 


8.5 ± 1, 


.7 


8.9 it 0.5 




5 6+0-2 
0.0_(, 2 


4.7 ± 1, 


.0 


5.6 it 0.3 


^-'-'-0.8 




1.4 ±0, 


.3 


2.9 it 0.3 


5 "S+1-0 


4 4+0-2 


2.6 ±0, 


.4 


2.1 it 0.1 


6 4+2-3 
°-*-1.8 


9 O+0.3 


2.2 ±0, 


.7 


1.6 it 0.2 



"Low signal-to-noise within the [0.15-l]r5oo aperture. T.y.soo is estimated 
using temperature within [0.15-0.5]r5oo and scaled using equation 8. 

'^SPT-CL J0551-5709 is coincident with Abell S0552. The gas temperature 
is estimated through equation 8. 



average ratio of the SZ- to the X-ray-mass is 0.89 ±0.06. 

Overall, we consider the agreement between the X-ray 
and SZ mass estimates reasonable given the ~ 15% sys- 
tematic uncertainty on the mass estimates quoted in VIO. 
However, there is some evidence that the SZ masses in 
VIO are biased low by ^--^10%. 

This difference could have several explanations. For 
example, the semi-analytic gas model used to calibrate 
the simulations in VIO could have a redshift evolution 
that differs from self-similar evolution. There could also 
be differences in the gas profiles at large radii that cause 
systematic differences between the measured Yx-based 
mass and the derived mass from the M-^ scaling in VIO 
that were not included in the simulations. It should also 
be noted that the AI^qq-Yx relation used here has been 
calibrated using X-ray hydrostatic mass estimates. 

Numerical simulations have suggested that hydrostatic 
mass estimates could be biased low by ~10% due to 
additional non-thermal pressure support of the gas in 
clusters (e.g., Lau et al. 2009; Meneghetti et al. 2010). 
However, comparisons of hydrostatic and weak lensing 
mass estimates of low redshift (z < 0.3) clusters show 
agreement at the <9% level and suggest no significant 
bias (e.g., Zhang et al. 2008; Vikhlinin et al. 2009a; 
Zhang et al. 2010). Additionally, the M^qq-Yx rela- 
tion assumes self-similar redshift scaling. The details of 
redshift-dependent departures from this scaling are not 
constrained observationally, although the results from 
simulations provide a 5% upper limit for any evolution 
of the amplitude of this relation between z = and 0.6. 

Regardless of the exact reasons, if the SZ mass esti- 
mates are biased low, there would be consequences for 
the cosmological constraints in VIO. It would gener- 
ally favor a lower normalization of the SZ significance- 
mass relation and hence a larger value of cts; however, 
these results should still be within the 68% confidence 
region in VIO due to the priors imposed on the parame- 



ters defining the SZ significance-mass relation which ef- 
fectively include the 15% systematic uncertainty for the 
masses. The cosmological implications of including the 
X-ray mass estimates with the VIO results will be the 
subject of a future paper. 

5. SZ DATA ANALYSIS 
5.1. Calculating Ysz 

The SPT observations, data reduction, and map- 
making used in this work are identical to VIO, and are 
outlined in §2.1. The analysis differs from VIO at the 
point where the SPT maps are spatially filtered to iden- 
tify the cluster candidates. In our case, the cluster candi- 
dates have already been identified, and we instead want 
to calculate the SZ-inferred integrated Comptonizatioii, 
Ygz, of each cluster. We calculate Ysz from the SPT 
150 GHz maps by spatially filtering them using a filter 
motivated by the X-ray measurements of each cluster. 

In VIO, cluster candidates are identified by spatially 
filtering the SPT maps with a matched filter (Haehnelt 
& Tegmark 1996; Melin et al. 2006). This is done in the 
Fourier domain by multiplying the map by 



B{K,ky)S{\k\ 



B(kxi ky)'^ Nastroi\k\) + Nnoise{kxj ky) 



(9) 



where "0 is the matched filter, B is the response of the 
SPT instrument after timestream processing to signals 
on the sky. S is the assumed source template, and the 
noise power has been broken into astrophysical (Nastro) 
and noise {Nnoise) components. The noise covariance 
matrix Nnoise includes contributions from atmospheric 
and instrumental noise, while Nastro includes power from 
primary and lensed CMB fluctuations, an unresolved SZ 
background, and unresolved point sources. In VIO, the 
source template, S, was constructed from a projected 
spherical /?-model. 
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TABLE 3 

Cluster masses from A^soo — Yx and A^sqo — Tx relations 



Name 


z 








A^500,sz,yio'' 






kpc 


10" M0 


10"Mq 


1O"M0 



JOOOO-5748 0.74 

J0509-5342 0.4626 

J0516-5430 0.2952 

J0528-5300 0.7648 

J0533-5005 0.8810 

J0546-5345 1.0665 

J0551-5709 0.4230 

J0559-5249 0.6112 

J2331-5051 0.5707 

J2332-5358 0.32 

J2337-5942 0.7814 

J2341-5119 0.9983 

J2342-5411 1.08 

J2355-5056 0.35 

J2359-5009 0.76 



950 


5.32 ± 


1062 


5.43 ± 


1463 


11.84 ± 


775 


2.97 ± 


656 


2.06 ± 


840 


5.33 ± 


948 


3.56 ± 


1043 


6.40 ± 


972 


4.70 ± 


1134 


5.66 ± 


1046 


7.43 ± 


847 


5.06 ± 


647 


2.47 ± 


1014 


4.18 ± 


816 


3.45 ± 



1.16 6.74 ±5.14 

0.60 6.71 ±3.40 

1.25 12.34 ±2.13 

0.89 3.05 ±4.00 

0.53 2.25 ±1.71 

0.62 5.25 ±1.13 

0.43 3.00 ±0.85 

0.54 7.07 ±1.49 

0.51 4.91 ± 1.36 

0.48 6.69 ±1.08 

0.90 8.10 ±2.18 

0.66 5.20 ±1.69 

0.32 2.39 ±0.63 

0.43 4.26 ±1.57 

0.67 4.97 ±2.61 



2.89 ±0.61 ±0.41 
4.26 ±0.74 ±0.60 
6.48 ±0.95 ±1.13 
2.83 ±0.60 ±0.38 
2.71 ±0.56 ±0.37 
3.25 ±0.51 ±0.44 
4.10 ±0.75 ±0.58 
5.03 ±0.74 ±0.70 
4.63 ±0.73 ±0.66 
5.19 ±0.85 ±0.83 
6.32 ± 0.84 ± 0.97 
4.05 ±0.58 ±0.63 
2.65 ±0.50 ±0.37 
4.17 ±0.80 ±0.63 
3.32 ±0.60 ±0.46 



Note. — Masses estimated from X-ray mass-proxy relations in equations 1 and 5. 
SZ derived masses from VIO arc shown for comparison. 

^Estimated using the Afsoo — Yx relation (equation 5). 

''Quoted masses from VIO include statistical and systematic uncertainties and 
have been scaled with our adopted h. 



In this work, we instead use a SZ source template mo- 
tivated from X-ray measurements of each cluster. The 
SZ brightness is proportional to the line-of-sight integral 
of electron pressure, or density times temperature. The 
profile is assumed to match the product of the best-fit gas 
density profile to the X-ray measurements of each clus- 
ter (equation 7), and the universal temperature profile 
of Vikhlinin et al. (2006) (equation 12). These profiles 
are multiplied together to give the radial pressure pro- 
file, and projected onto the sky by doing a line-of-sight 
integral through the cluster. The radial pressure profile 
is truncated at 3 x r5oo, where rsoo is estimated using 
equation 5 and given in Table 3. In constructing the 
spatial filter, we only need to know the spatial shape of 
the source, so the SZ model is normalized to unity. 

For each of the fifteen clusters, we construct a spatial 
filter using equation 9 and a source model, as described 
above. The spatially filtered SPT maps are a measure 
of the normalization of each source model at the cluster 
location. Using the SPT maps alone to determine the 
cluster location would bias the Ysz measurements high 
because of the freedom to maximize the SPT significance 
by position. Therefore we use the X-ray measured posi- 
tion as a prior on the cluster location. We define the un- 
certainty in the SPT normalization of the source model 
as the standard deviation of the spatially filtered map 
within a 90' band in declination around the cluster. The 
SPT maps are calibrated in units of Kca/b, the equiv- 
alent CMB temperature fluctuation required to produce 
the same power fluctuation. The SPT normalization is 
converted to Comptonization using the relation 

AT^yTcMB(^x^;^-4j[l + 6{x,T,)], (10) 

where AT is the measured decrement in units of Kcmb , 
y is the Comptonization, Tqmb is the CMB blackbody 
temperature of 2.725 K, 2: = hv /kTcMBi a-nd 5{x^Te) 
accounts for relativistic corrections to the SZ spectrum 
(Itoh ct al. 1998; Nozawa et al. 2000). For the frequency 



dependent terms in equation 10, we calculate their SPT 
band averaged value, which would effectively correspond 
to an observing frequency, of 152.9 GHz. The aver- 
age SPT band is measured from Fourier Transform Spec- 
troscopy measurements of a sample including more than 
90% of the SPT 150 GHz detectors. 

For the relativistic correction factor, 5{x,Tf.), we as- 
sume an electron temperature of Tx,r50o/l-lli where 
Tx,r500 is given in Table 2 and the factor of 1.11 is the av- 
erage ratio of the X-ray spectroscopic and mass- weighted 
temperature measured in Vikhlinin et al. (2006). The 
relativistic correction factor is only a weak function of 
temperature, between 3 to 10 keV it varies by ~ 4%, 
and we expect this assumption to negligibly affect our 
results. 

The integrated Comptonization, Ysz,cyh is calculated 
for each cluster by integrating its source model over solid 
angle, Ysz,cyi oc / y{6)dil, normalized to the best-fit SPT 
Comptonization, y. To more easily compare to the X-ray 
measurements, we convert our measurements to units of 
M0keV, and define Ysz,cyi as. 

Ysz,cyi = (^t^^H^^^^ D\ J y{9)dn, (11) 

where Da is the angular distance to the source, (Jt is 
the Thomson cross-section, nie is the electron mass, rup 
is the proton mass, c is the speed of light, fie is the mean 
molecular weight per free electron. In Table 4, we give 
ysz,cyi integrated out to the angular radius correspond- 
ing to r5oo as determined from the X-ray measurements 
and given in Table 3. The uncertainty in Ygz is calcu- 
lated as the quadrature sum of the uncertainty in the 
SPT calibration, the measured SPT normalization of the 
source model, and an additional uncertainty in the as- 
sumed source model, which will be discussed in §5.3. 

5.2. Spherical Deprojection 

Both the X-ray and SZ observations are measuring a 
projected signal that is proportional to the integrated 
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TABLE 4 

Measured cylindrical and spherical Ysz 



Name 


z 


ySZ.cyl.500 
10" Mq keV 


ySZ.sph.bOO 

1014 Mq keV 


JUUUU-0 / 4(5 


0.74 


3.0 it 0.6 


2.4 it 0.5 


JUoUy-oo42 


0.4626 


3.9 it 0.8 


3.1 it 0.6 


JUolO-04ciU 


0.2952 


15.3 it 2.4 


11.6 it 1.8 


JUozo-OoUU 


n vfi/i Q 
U. ( 04o 




o 1 -1- n /I 
z. i it U.4 


J0533-5005 


0.8810 


1.8 ± 0.5 


1.5 it 0.4 


J0546-5345 


1.0665 


4.6 it 0.8 


3.5 it 0.6 


J0551-5709 


0.4230 


3.7 it 0.8 


2.6 it 0.6 


J0559-5249 


0.6112 


7.6 ± 1.1 


5.8 it 0.9 


J2331-5051 


0.5707 


5.0 ±0.8 


3.6 it 0.6 


J2332-5358 


0.32 


5.6 ± 1.0 


4.9 it 0.9 


J2337-5942 


0.7814 


8.3 it 1.0 


7.1 it 0.9 


J2341-5119 


0.9983 


4.9 it 0.7 


4.1 it 0.6 


J2342-5411 


1.08 


2.3 it 0.5 


1.9 it 0.4 


J2355-5056 


0.35 


2.6 it 0.6 


2.2 it 0.5 


J2359-5009 


0.76 


2.1 ± 0.6 


1.7it0.5 



gas properties in a line of sight through the cluster. This 
projection has different physical dependencies between 
X-ray and SZ observations that must be considered. To 
ease this comparison we deproject each measurement so 
that they correspond to a physical observable that is pro- 
portional to a spherical volume integral of each signal. 

The Yx estimates are deprojected as described in §3.4. 
For X-rays, the effect of the deprojection is smaller than 
for the SZ, because the X-ray flux is proportional to 
and only weakly dependent on temperature. This effect 
decreases the contribution to the X-ray signal from large 
radii where the density is lower. Also, Yx will be pro- 
portional to a X-ray spectroscopic weighted temperature 
from gas between 0.15r5oo and rsoo, while Ysz is related 
to the gas mass weighted temperature. We do not for- 
mally account for this difference in calculating either Yx 
or Ysz- However later in this section we will discuss this 
effect in interpreting their comparison. 

We calculate Ysz,sph by spherically deprojecting the 
Ysz,cyi measurements in Table 4. The SZ signal is pro- 
portional to the line of sight integral of the electron den- 
sity, TT-e, times temperature, Tg. For each cluster, we as- 
sume the density profile derived from the X-ray imaging 
analysis and assume a temperature profile of a universal 
form. 



T{r) ^ To 



(a;/0.045) 



1.9 



0.45 



1 



(12) 



(a;/0.045)i-9 + 1 (1 + (.t/0.6)2)0-45 ' 

where x = r/r^Qo (see Vikhlinin et al. 2006). We then 
define Ysz,sph as: 



Y, 



SZ,sph 



SZ.cyl 

c 



Y. 



IspHnr)ne{r)dV 



SZ,cyl 



Icyinr)ne{r)dV ' 



(13) 



where C is the ratio of the integrals of pressure in a 
cylindrical and spherical volume through the cluster. For 
these integrals we use an integration radius of rsoo and 
truncate the density and temperature radial models at 
Sj'soo- We note that varying the truncation radius from 
between 3-5r5oo changes our measurement of Ysz in Ta- 
ble 4 by less than 1%. For the r^oo aperture, the me- 
dian and standard deviation of C across the sample is 
1.23 ± 0.08. This is consistent with the value of 1.203 
that is calculated assuming the universal pressure profile 
from Arnaud et al. (2009). 



As noted earlier in this section, even after the above 
spherical deprojection, Yx and Ysz.sph are not directly 
comparable because of their different weighting of the 
electron temperature. Yx is proportional to the X-ray 
spectroscopic weighted temperature, Tx, while Ysz.sph 
is proportional to the mass weighted temperature, Tmg- 
The size of this difference has been estimated by several 
authors from X-ray measurements. In Vikhlinin et al. 
(2006), they estimate Tx/T^g = l-H for a sample of 
relaxed massive clusters with high angular resolution X- 
ray temperature measurements between 70 kpc and rsoo • 
This would imply that Ysz,sph/Yx = 1/1.11 « 0.90. In 
Arnaud et al. (2009), a similar analysis was performed 
for a sample that included both relaxed and un-relaxed 
clusters, and it was estimated that Ysz.sph/Yx = 0.924± 
0.004. For the measured density profiles and assumed 
temperature profiles for our sample, we would expect an 
average ratio of Fg^^sp/j/Yx = 0.91 ± 0.01. To compute 
this expected ratio we have used the "spectroscopic-like" 
temperature (e.g., Mazzotta et al. 2004), given the den- 
sity and temperature profiles used here. 

In §6.2 we compare the above ratios, predicted purely 
from X-ray observations, to the ratio including the 
Ysz,sph calculated from the SZ data as described in §5.1 
and deprojected as described in this section. 

5.3. Model Uncertainty 

The integrated Comptonization, Ysz, inferred from the 
SZ data depends on the assumed spatial model through 
the construction of the spatial filter, and the volume inte- 
gral through the deprojection the SZ data. For a typical 
cluster, the X-ray data constrains the cluster density pro- 
file with high signal-to-noise out to ~ ''soo/S with no in- 
formation on the temperature profile. Since a significant 
amount of the SZ signal is coming from larger radii than 
this, we would like to estimate how much uncertainty 
our assumed profile is adding to the Ysz estimates. To 
help do this, we calculate the Ysz of each cluster assum- 
ing the universal pressure profile of Arnaud et al. (2009), 
variations of the temperature profile in equation 12, and 
an isothermal model. 

Arnaud et al. (2009) measure a universal pressure pro- 
file from X-ray measurements for a representative set of 
local massive clusters. These measurements were deep 
enough to constrain the cluster density and temperature 
profiles out to rsoo in each cluster. They find that their 
sample is well fit by a universal pressure profile that is 
defined only by M500. For each cluster in our sample, we 
re-calculate Ysz,sph assuming the universal pressure pro- 
file from Arnaud et al. (2009) and the X-ray measured 
values for r^^o and Afsoo given in Tables 4 and 3, respec- 
tively. Comparing these Ysz.sph estimates to the values 
given in Table 4, we find that this ratio has a mean and 
standard deviation of 1.01 ± 0.15 averaged over the sam- 
ple. There is one significant (> 20%) outher, SPT-CL 
J0516-5430, whose Ysz.sph deviates by a factor of 0.59. 
This cluster is a major merger, and the observed cluster 
profile is noticeably poorly fit by the Arnaud et al. (2009) 
pressure profile, which does not capture the disturbed 
distribution of the central gas in this cluster. Leaving 
this cluster out of our comparison, the average ratio be- 
comes 1.03 ± 0.10. Therefore we see no detectable bias 
when assuming the Arnaud et al. (2009) pressure profile, 
but these results suggest that there could be an addi- 
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tional ^ 10% uncertainty in our Ysz.sph measurements 
from our assumed pressure profile. 

As an additional test, we vary the outer slope of the 
Vikhlinin et al. (2006) temperature profile in equation 12, 
letting the exponent on the term in the right hand side 
of the denominator vary between 0.16 to 0.80, from its 
starting value of 0.45. This range of values matches the 
full range of effective slopes of the temperature profile at 
r5oo as measured in the sample of Vikhlinin et al. (2006). 
Calculating new Ysz.sph estimates assuming these new 
temperature profiles, we find they change the Ysz.sph 
values given in Table 4 on average by a factor of 1.09 ± 
0.06 and 0.91 ± 0.05 for the exponent values of 0.16 and 
0.80, respectively. While we consider such a significant 
systematic shift in the temperature profile unlikely across 
the whole sample, this level of variation in temperature 
profiles could contribute added statistical uncertainty in 
the Ysz.sph estimates. 

As a limiting case, we also recalculate Ysz,sph assuming 
that the cluster is isothermal and with a density profile 
corresponding to the best-fit profile to the X-ray mea- 
surements. Comparing Ysz,sph for the isothermal profile 
to the values given in Table 4, we find that this ratio 
has a mean and standard deviation of 1.12 ± 0.08 aver- 
aged over the sample. This is certainly an extreme case 
because of the abundance of evidence for the temper- 
ature profile in clusters decreasing significantly by rsoo 
(e.g., Vikhlinin et al. 2006; Arnaud et al. 2009). Recent 
SuzaJcu-observations of nearby clusters at large radii also 
indicate a drop of gas temperature by a factor 3 at 
'"200 (e.g., Fujita et al. 2008; George et al. 2009; Bautz et 
al. 2009; Reiprich et al. 2009; Hoshino et al. 2010). We 
also note that the temperature profile with an exponent 
of 0.80 assumed in the previous paragraph would have a 
gas temperature that dropped by a factor of 5 at r2oo, 
so even this would seem like an extreme case. 

Overall, we are encouraged that the variation in 
Ysz,sph is found to be ^ 10% when assuming a broad 
range of different temperature and pressure profiles found 
in the works of Vikhlinin et al. (2006) and Arnaud et 
al. (2009). Therefore we conservatively assume an ad- 
ditional 10% statistical uncertainty in our Ysz.sph esti- 
mates from any assumed model uncertainty, which we 
have added in quadrature to the uncertainties given in 
Table 4. 

5.4. Point Source Contamination 

Astrophysical point sources in the direction of the clus- 
ter can potentially fill in the SZ cluster decrement and 
bias the integrated Comptonization low. There are gen- 
erally two populations of sources that can contaminate 
the SZ signal: dusty or radio bright sources. In VIO, the 
potential contamination from both were discussed, and 
neither is expected to significantly bias the SZ measure- 
ments at 150 GHz averaged over the sample. We review 
some of those conclusions here. 

In the Appendix, we discuss radio detections at 843 
MHz by the Sydney University Molonglo Sky Survey 
(SUMSS) towards the clusters in this work. The majority 
of radio sources in clusters have been found to have steep 
spectra with a < —0.5 (where S (x ly") (e.g.. Coble et al. 
2007; Lin et al. 2009). For example. Coble et al. (2007) 
find a median spectral index of —0.72 between 1.4 and 
28.5 GHz for radio sources towards a sample of massive 



clusters ranging from 0.14 < z < 1.0. In VIO, they noted 
that a typical cluster would suffer a decrease of = 1 
for a 2 mJy (5 mJy) source at 150 GHz. located at 0.5' 
(1') from the cluster center. Assuming a spectral index 
of —0.72, a 2 mJy (5 mJy) source at 150 GHz would be 
~83 (210) mJy at the SUMSS observing frequency of 843 
MHz. As detailed in the Appendix, no source has been 
detected above either threshold within 1' of any cluster 
in this work. However, for any individual cluster radio 
source, extrapolating its flux from radio frequencies to 
the SPT observing frequency of 152.9 GHz is difficult 
because of the broad range and frequency dependence of 
the spectral indices that is typical for these sources (e.g.. 
Coble et al. 2007; Lin et al. 2009). In practice, extrapo- 
lating the fiux of any SUMSS detected source to the SPT 
observing frequency will have factors of a few uncertainty 
without further observations spanning an intermediate 
range of frequencies. With this caveat in mind, there ap- 
pears to be no source strong enough to significantly bias 
the SZ fiux estimate for any cluster in this work. 

Including every source detected within 0.5' of each 
cluster center in our sample, the cumulative SUMSS flux 
is 224.5 mJy. If we assume a spectral index of -0.72, this 
would imply an average flux of ^0.35 mJy per cluster 
at 152.9 GHz. This would correspond to an average de- 
crease in the SZ flux of ^2%, assuming that these sources 
represent an overdensity to the background population. 

The average radio source contamination in SZ surveys 
has been recently estimated in simulations of the mi- 
crowave sky by Sehgal et al. (2010). These simulations 
were motivated by observations by Lin et al. (2009) mea- 
suring the radio source population characteristics in low 
redshift {z < 0.25) X-ray detected clusters. For the mass 
limit of this work, M500 w 3 x IO^^H^^Mq, these stud- 
ies predict that in <1% of clusters there would be radio 
source contamination large enough to effect the SZ flux 
measurement at the > 20% level, with this result largely 
independent of redshift. VIO estimates a similar rate of 
contamination using a radio source count model (de Zotti 
et al. 2005) and the measured over-abundance of radio 
sources near clusters (Coble et al. 2007). Overall, the 
combination of the above results lead us to not expect 
any significant radio source contamination. 

Emission from sub-millimeter bright galaxies can also 
contaminate the SZ signal. These sources are typi- 
cally dusty star forming galaxies that are very lumi- 
nous but highly obscured, such that their luminosity 
is peaked at infrared (IR) wavelengths. The number 
counts and fluxes of these sources have been character- 
ized between 0.5 - 2.0 mm wavelength by several experi- 
ments (e.g., Coppin et al. 2006; Devlin et al. 2009; Vieira 
et al. 2009; Austermann et al. 2010), and these mea- 
surements reasonably match simple analytic models de- 
scribing their source population distribution (Lima et al. 
2010b). These sources can largely be approximated as a 
random background that contributes additional signal in 
the SPT maps, and we have explicitly accounted for them 
in our matched filter in equation 9. However, this imple- 
mentation did not account for any emission that could 
be correlated with clusters, or additional noise from the 
background objects being gravitationally lensed by the 
cluster. The former was argued to be insignificant in VIO 
because: the IR overdensity towards clusters is expected 
to be ^ 1 mJy at 150 GHz even out to highest redshift 
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clusters in our sample, and observational measurements 
of the 100 /im flux towards local clusters {z ^ 0.2) are 
too low to significantly bias the SZ flux measurements 
even allowing for a strong redshift evolution of these 
sources. The latter has been characterized by Lima et 
al. (2010a) to increase the flux noise towards clusters by 
~60%, which is at a level such that it would be <3% of 
the SZ flux for the mass range of clusters in our sam- 
ple. We note that because lensing is a surface brightness 
conserving process, this latter effect would not bias our 
results when averaged over a large sample of clusters. 

One noteworthy exception where point source contami- 
nation is evident is the cluster SPT-CL J2332-5358. This 
cluster was the only candidate in VIO that was also co- 
incident with a point source detected in the SPT 220 
GHz maps at > 5a. Similar 220 GHz detected point 
sources have a source flux, S, that scales with frequency 
as 5 oc J/" with a = 3.3 ± 0.7 between 150 and 220 GHz 
(Vieira et al. 2009). Extrapolating its measured 220 GHz 
source flux would predict a flux of 11.0 ± 2.7 mJy at 150 
GHz. To the 150 GHz SPT maps, we subtract a point 
source of this brightness, after convolving it with the 150 
GHz SPT transfer function. We then repeat the same 
cluster extraction method as in VIO to calculate a new 
cluster position and SPT-significance, ^, which is given 
in Table 1. This point source corrected 150 GHz map 
is used when calculating the Ysz given in Table 4, and 
the added uncertainty in Ysz from the uncertain point 
source spectral index is added in quadrature with the 
other uncertainties outlined in §5.1. 

6. SCALING RELATIONS 

Here we discuss the X-ray and SZ scaling relations for 
the sample. We consider the integrated SZ flux Ysz 
for both cylindrical and spherical volumes of integration. 
The integrated cylindrical estimate is deprojected to a 
spherical estimate as described in §5.2. 

To measure the slope and normalization of the scaling 
relations we perform linear regression in the presence of 
intrinsic scatter. We use the Bayesian method of Kelly 
(2007) and maximize the probability of a linear model 
in log-log space, accounting for measurement uncertain- 
ties in both axes. For the Ysz scaling relations below 
we fit the relations with another method, accounting for 
the ^ selection. We investigate this fitting method using 
simulated cluster samples as detailed in §6.2. We have 
assumed self-similar E{z) scaling for these scaling rela- 
tions unless otherwise stated. The best fit parameters 
for these relations are listed in Table 5. 

6.1. Mg - Tx, Lx - Tx and Lx-Msoo 

The scaling between the ICM gas mass {Mg) and the 
spectroscopic temperature (Tx) is an important test of 
the properties of this SZ-selected sample. We want to 
know if the SZ selection biases the observables in any way 
or if they are consistent with those from X-ray selected 
samples. The values of Tx and Mg within r^oo used here 
were derived as described in §3.3 and §3.4 and are given 
in Table 2. This relation can provide clues about the 
high-z behavior of the M-Tx and M-Mg scaling rela- 
tions (equations 1 and 4). To accurately constrain the 
evolution of these relations, however, either hydrostatic 
masses or masses obtained via an independent measure- 
ment are needed for high-z clusters. 



We fit the relation 

to the data as described above, accounting for intrinsic 
scatter. Figure 1 shows the Mg — Tx relation with the 
best fit power law (solid line). 

To compare our results with previous studies using X- 
ray selected samples we use the gas mass and temper- 
ature data for the low-z (median z « 0.05) and high- 
z (median z w 0.48) as measured by Vikhlinin et al. 
(2009a) . We choose these samples since the data analysis 
is analogous to ours. We fit the scaling relation to these 
samples using the same method. For our SPT-selected 
sample (median z = 0.74) we find A = 13.65 ± 0.10 and 
B = 1.95 ± 0.66 (see Table 5). When fitting the low-z 
sample, we find A = 13.54 ± 0.02 and B = 1.66 ± 0.08, 
at 5 keV. For the high-z sample from Vikhlinin et al. 
(2009a), the best fit parameters are A = 13.66 ± 0.03 
and B = 1.64 ± 0.21 in even better agreement with our 
findings. This offset between low-z and high-z samples 
indicate deviations from self-similarity. The slope ap- 
pears steep in our relation, although consistent with both 
the low- and high-2 X-ray selected samples as well as the 
self-similar slope, B = 1.5. A positive trend of the gas 
mass fraction with mass has been observed previously 
(e.g., Mohr et al. 1999; Vikhhnin et al. 2006) and is likely 
what is causing the steeper than self- similar slope of the 
relation. 

Numerical simulations (e.g. Kravtsov et al. 2005) also 
show an increasing gas mass fraction with redshift, pre- 
sumably due to a different distribution of mass progeni- 
tors at high redshift. The SPT-selected cluster sample is 
consistent with this scenario. The normalization agrees 
well with that of the high-z sample of Vikhlinin et al. 
(2009a) while still consistent with the low-z sample at 
the 1(7 level. We cannot rule out that this behavior is 
partly related to deviations from self-similarity in the 
M-Tx relation at high redshift. However, the results 
are consistent with previous findings from X-ray selected 
cluster samples. 

We also investigate the luminosity-temperature rela- 
tion (e.g.. Edge & Stewart 1991; Markevitch 1998) for 
the sample where we compare the luminosity in the 0.5- 
2.0 kcV band within rsoo, including the core, Lx, to the 
temperature estimate within ^500, Tx, with the core ex- 
cised. We fit the relation 

i^(.)-L.=.10-(^)', (15) 

and find a shallow slope, B ~ 1.92 ± 0.60 (see Table 5), 
marginally inconsistent with most previous works which 
typically find a slope close to 3 (e.g. Mantz et al. 2009). 
For completeness, we also estimate the bolometric lumi- 
nosity by extrapolating the spectrum for all frequencies 
using the core excised temperature. This introduces ad- 
ditional uncertainties on the luminosity due to the un- 
certainty in the temperature estimate. For the Lboi-Tx 
relation, we find a steeper slope of B = 2.31 ± 0.85, con- 
sistent with previous work. The best fit relation is shown 
in Figure 2 (left). 

We investigate the luminosity- mass relation (e.g., 
Reiprich & Bohringer 2002) and adopt the best-fit E{z)- 
scaling found by Vikhlinin et al. (2009a) to fit for the 
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Fig. 1. — Mg — Tx relation with self-similar _E{^)-scaling. The data from the SPT sample are shown as circles with error bars and the 
solid line shows the best fit scaling relation. For comparison we also show the best fit relations using the low-z (dashed line) and high-2 
(dotted line) samples from Vikhlinin et al. (2009a). These data were fitted using the same method described here. SPT-CL J0551-5709 is 
shown as a red square and SPT-CL J0516-5430 is shown with a blue triangle. 



relation 

Lx = 10-(^^)'. (16) 

The best fit relation is shown in Figure 2 (sohd hne) com- 
pared to the best fit in Vikhhnin et ah (2009a) (dashed 
fine). We find a slope of _B = 1.16 ± 0.20, shallow com- 
pared to the B ~ 1.61 ± 0.14 found by Vikhlinin et al. 
(2009a), while the normalization agrees very well near 
the median mass of the sample at 5 x 10^^ Mq. From 
the relation in Vikhlinin et al. (2009a) we expect a nor- 
malization oi A = 44.25 at 5 X lO^'' Mq, which is in good 
agreement with our measured value oi A = 44.28 ± 0.07. 
The slope of the fit is mostly driven by SPT-CL 0516- 
5430 (blue triangle) which has a high derived mass for 
its measured luminosity, possibly caused by the observed 
merging activity in the cluster. Excluding this object 
from the fit changes the slope to B = 1.45 ± 0.29, con- 
sistent with Vikhlinin et al. (2009a). 

In summary, we find a good agreement between the 
X-ray observables measured for this SZ-selected sample 
when comparing our results to X-ray selected samples 
where similar analysis methods were used. 

6.2. Ysz scaling relations 
6.2.1. Fitting method 

The selection of clusters above a fixed S/N threshold 
can bias the SZ-X-ray scaling relations. The SZ flux, 
is correlated with the SZ S/N on which the sample 
is selected and ignoring this effect will lead to bias in the 
Yg^-scaling relation parameters. The steepness of the 
cluster mass function makes this bias more of a problem 
since the number of clusters in a given mass bin grows 
exponentially towards lower masses. This causes more 
clusters to be near the selection threshold where the ef- 
fect is most prominent. The measurement uncertainty in 
Ysz and the intrinsic scatter in the scaling relation cause 
clusters with low mass (or Yx) to scatter over the selec- 
tion threshold. This is visible as a tail towards low mass 
(or Yx) near the Fgz-selection threshold in the scaling 
relation plots. If ignored, this bias will lead to a biased 
slope and a higher normalization near the threshold in 
the fitted scaling relation parameters. Below we describe 



our method to account for the SZ selection in our scaling 
relation fits. 

We estimate how the selection cut in ^ translates to 
a selection in Ysz,sph using simulated SPT observations. 
We repeat the procedure used in VIO to estimate the 
SZ selection function and SZ-significance mass relation, 
where observations of simulated SZ maps are added to 
the dominant sources of astrophysical confusion and in- 
strumental noise, mock observed, and processed through 
the SPT analysis pipeline. In contrast from VIO, we keep 
track of the predicted value of Ysz,sph for each cluster, 
which is estimated using the Ygp/i — -^500 scaling rela- 
tion from Arnaud et al. (2009) and the simulated clus- 
ter's mass, AfsoQ. Also, we measure ^ for each simulated 
cluster in the same way as we do for the SPT observa- 
tions, where we record the maximum significance S/N 
over different spatial filters and in a single realization of 
the astrophysical confusion and noise. 

With these results we find a best-fit relation 

( 2xl0iS/ekev ) 
with intrinsic scatter in logYgz of 0.10. Comparing this 
result with the measured ^ and Ysz for our sample we 
find that the simulation-based relation is offset slightly 
high in ^. When fitting the normalization of the above 
relation with the observed data we find a best fit value of 
5.56 ±0.31, about 6% lower than the simulations suggest. 
The reason for this offset may be related to systematic 
uncertainties in the simulation of SPT observations. We 
adopt this best fit normalization from the observed data, 
and use the above relation to construct a selection func- 
tion in logYsz- We estimate the selection probability as 
an error function in \ogYsz at this threshold with the 
width set by the intrinsic scatter in this relation. To 
account for the uncertainty in the simulations we use in- 
clude a 10% gaussian uncertainty on the normalization 
and a 20% gaussian uncertainty on the intrinsic scatter, 
which we marginalize over for our results. 

The Ysz scaling relation fit is performed as follows. 
First, we calculate the probability of measuring a certain 
Ysz for a given M or Yx and a given scaling relation 
with log-normal scatter. This is then convolved with the 
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kTlkeV] M,^_y^lM,^J 

Fig. 2. — L),ai — Tx (left) and Lx — Afsoo.yx (right) relations for the sample. For the Lf^^i — Tx relation we have applied self-similar 
E{z) scaling. For the Lx — A^soo.Vx relation we have applied i?(^)^ *^-scaling to enable comparison with Vikhlinin et al. (2009a) shown 
with a dashed line. Solid lines represent the best fit relation. SPT-CL J0551-5709 is shown as a red square and SPT-CL J0516-5430 is 
shown with a blue triangle. 



measurement uncertainty in Ysz- The Ysz selection cut 
is applied by multiplying the probability of measuring 
Ysz with an error function in logYsz as described above. 
The likelihood is then calculated as the product of prob- 
abilities for all clusters and is maximized to obtain the 
scaling relation parameters. 

The uncertainties on the calculated Ysz are correlated 
with the measurement error in Yx and M^qq^Yx ■ The 
r5oo radius used for integration of Ysz is estimated from 
M5oo,Yx and any measurement scatter in Yx will con- 
tribute to scatter in Ysz of similar magnitude. This will 
lead to scatter along the scaling relation and for this rea- 
son we ignore these effects. We fit for the intrinsic scatter 
in the Xsz-Afsoo.yx and Ysz-Yx relations while noting 
that the scatter in Ysz and Yx with the true cluster mass 
will also be partially correlated. 

To test our fitting method for potential biases we run 
it on mock samples drawn from a fiducial mass function. 
We perform this test both in the limit of a sample with 
1000 clusters and for many samples of 15 clusters. To es- 
timate a, Ysz for the mock samples we assume a Ysz-M 
relation with log-normal scatter and convolve this with 
a linear measurement scatter of 0.4 x 10^^ Mq. We as- 
sume that we know the relation between ^ and Ysz and 
its scatter to reproduce the sample selection function. In 
the limit of a large sample of clusters, we find no measur- 
able bias in either the fitted slope or normalization. For 
1000 generated mock samples of 15 clusters we find that 
the scatter in normalization and slope between the sam- 
ples and the average measurement error on the parame- 
ters over the samples are consistent. The normalization 
is reproduced to within 0.1% of the input value while 
we find that the slope is biased low by 2.0%. The aver- 
age measurement error on the slope from the 15 cluster 
samples is 13.5% so this bias is of little significance. 

6.2.2. Ysz-Yx relation 

The relationship between Ysz and Yx is determined 
by details of the gas and temperature distribution in the 
cluster. It is effectively measuring a relationship between 
the mass weighted and X-ray spectroscopic weighted 
temperature. In Arnaud et al. (2009), X-ray measure- 
ments of both relaxed and unrelaxed clusters were used 



to approximate a Ysz-Yx relation, for a spherical inte- 
gration to r.goo- Their work found a relationship with 
a normalization of 0.924 ± 0.004, implying a lower mass 
weighted temperature. This result is consistent with pre- 
vious X-ray measurements (e.g., Vikhlinin et al. 2006), 
and our expectation of 0.91 ± 0.01 for the measured gas 
density profiles and the assumed temperature profile in 
this work. Hydrodynamical simulations predict a simi- 
lar ratio, but also find that X-ray measurements could 
overestimate the integrated pressure at some level (e.g., 
Nagai et al. 2007). This would imply that the X-ray 
predictions of the normalization of the Ysz-Yx relation 
could be biased high; for example, the results of Nagai 
et al. (2007) suggest a bias at the few percent level. 

Recently there has been some evidence that the SZ 
signal from clusters could be lower than this expectation 
from comparison of WMAP SZ observations of X-ray se- 
lected clusters (Komatsu et al. 2010), however similar 
analyses have led to different conclusions (Melin et al. 
2010). By fitting a normalization of the Ysz-Yx relation 
to the clusters in this work, we can test for a similar off- 
set, which could be indicative of differences in the real 
gas profiles from those we've assumed or some other sys- 
tematic bias in either the SZ or X-ray measurement. 

Using the method described in §6.2.1, we fit a scaling 
relation between Ysz and Yx of the form 

In Figure 3, we show the Ysz-Yx relation for both cylin- 
drical (left) and spherical (right) Ygzi denoted by Ysz,cyi 
and Ysz,sph respectively. In Table 5, we give fits to equa- 
tion 18 given several different assumptions, which are 
described and discussed further below. 

For both cylindrical and spherical Ysz we find a slope 
consistent with unity. The dashed line in the plots is the 
Ysz = Yx relation and the dotted line (right) shows the 
expected Ysz,sph/Yx ratio of 0.924 found in the work of 
Arnaud et al. (2009). We also fit the scaling relations 
with the slope fixed to 1 (see Table 5). The normal- 
ization of our Ysz.sph-Yx fit implies an average ratio of 
0.82±0.07, consistent with the prediction of Arnaud et al. 
(2009) at the 1.5(t level. The most significant outlier from 
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the Ysz,sph-Yx relation, SPT-CL JOOOO-5748, also has 
the lowest Ysz,sph/Yx ratio in our sample. If we exclude 
this cluster from the fit, we find an even better agreement 
with a normalization of 0.87 ±0.06. The marginally lower 
normalization found in this relation when compared to 
the X-ray prediction is expected at some level from hy- 
drodynamical simulations which predict that Ysz could 
be biased high when estimated using X-ray data. 

As an additional test, we fit the Ygz-Yx relation in- 
cluding a factor of £'(2;)'^, and fit for C while keeping the 
slope fixed to 1. We find C = —0.36 ± 0.58, consistent 
with no evolution. Further investigating this, we divide 
the sample in two redshift bins. The low-z bin consists 
of the 7 clusters at 2; < 0.7 while the high-z bin consists 
of the 8 clusters at z > 0.7. Fixing the slope to 1, wc find 
a ratio of 0.88 ± 0.12 for the low-z sample, and a ratio 
of of 0.72 ± 0.14 for the high-z sample. If we exclude the 
most significant outlier from the high-z sample, SPT-CL 
JOOOO-5748, we find a ratio of 0.83 ± 0.09. 

As discussed at the beginning of this section, the 
Ysz,sph/Yx ratio is expected to be lower than 1 due to 
the different weighting of the temperature in the two es- 
timates. As explained in §5.3, our measurement of Ysz 
is sensitive to the assumed temperature profile, whereas 
the Yx measurement is not. This causes the measured 
Ysz,sph/Yx ratio to change under different assumptions 
about the temperature profile. However, the expected 
ratio of the gas-mass weighted and the spectroscopic 
weighted temperatures changes by a similar amount as 
the temperature profile changes. The overall effect makes 
the measured ratio to be fairly independent to our as- 
sumptions of the shape of the temperature profile. This 
is particularly true for variations in the slope of the tem- 
perature profile at large radii that we tested in §5.3. 

Overall wc find a normalization consistent with our 
expectations given the gas profiles in this work, and be- 
lieve this result to be largely independent of our assumed 
temperature profile. Wc find no significant evidence for a 
redshift dependent evolution in the normalization of the 
Ysz-Yx relation, but note that the high-z sub-sample 
does favor a '^Icr lower normalization. This could be 
marginal evidence for redshift evolution, however further 
X-ray and SZ observations will be needed to make any 
statistically significant statements. 

6.2.3. Ysz-M5oo,Yx relation 

Finally, we investigate the relation between the SZ fiux, 
Ysz, and the X-ray derived mass, M^qq.Yx- This is not 
an independent result from the previous section because 
the X-ray masses are calculated directly from the Yx 
measurements and the Mr^QQ Yx relation in Vikhlinin 
et al. (2009a). However, it is useful in understanding the 
SZ mass calibration, and can be compared to previous 
measurements of this relation. 

Using the method of §6.2.1, wc fit a scaling relation 
between Ysz and A/5oo,Yx of form 

Ysz = 10-^ ( ) E{zfl\ (19) 

V3 X 10" Mq) ^ ' ^ ' 

In Figure 4 we show the Ysz — M^qq yx relation, and in 
Table 5 wc give the best-fit parameters to equation 19. 
The slope of this relation is found to be 1.67 ± 0.29 for 
Ysz,sph-M^QQ^Yx with intrinsic scatter of 0.09 ±0.05 (see 



Table 5). This slope is consistent with the self-similar ex- 
pectation of 5/3 and previous measurements (e.g., Bona- 
mente et al. 2008; Arnaud et al. 2009). In Arnaud et 
al. (2009), a slope of 1.790 ± 0.015 is measured for their 
predicted Ysz.sph using X-ray observables. 

Wc also fit the relation keeping the slope fixed at the 
value expected from X-ray studies, 1.79, and note that 
the normalization of the Ysz-M relation is lower than 
that of Arnaud ct al. (2009). The Arnaud et al. (2009) 
results imply a normalization A = 14.115 ± 0.003, using 
our adopted h while we find A = 14.03 ± 0.04, around 
a 2a offset. The best fit relation is shown in Figure 4 
(solid line) with the Arnaud et al. (2009) relation shown 
(dashed line) for comparison. 

It should be noted that part of the offset is due to 
differences in the mass estimates in our work and in 
the work of Arnaud et al. (2009). The normalization in 
the M — Yx relation (equation 5), used here to estimate 
cluster masses, is different from the normalization of the 
M — Yx relation used in Arnaud et al. (2009) to derive 
the Ysz-M5oo,Yx relation in that work. For our adopted 
h, wc find a mass, 

M5oo,Yx 4.83 ± 0.17 X lO" Mq, at 
Yx = 3x 10" Mq kcV from equation 5. If we instead 
use the M — Yx relation from Arnaud et al. (2009), we 
find Ahoo^Yx = 4.64 ± 0.12 x lO" Mq, at the same Yx- 
To account for this scaling relation offset we shift the 
normalization of the Arnaud et al. (2009) Ysz-Msqq^Yx 
relation down by a factor (4.83/4.64)1™ = 1.074, result- 
ing in ^ = 14.084 ± 0.003. This is reduces the offset to 
our measurement, A = 14.03 ± 0.04, to a similar level to 
what we found for the Ysz-Yx relation in the previous 
section. 

It is important to note that the intrinsic scatter in this 
scaling relation does not directly reflect the low scat- 
ter relationship between Ysz and the gravitational mass. 
The mass is derived directly from Yx and therefore in- 
cludes the scatter in the Yx-mass relationship. Inde- 
pendent measurements of the gravitational mass, e.g., 
through weak lensing or hydrostatic masses from X-ray 
data, are necessary to diagnose scatter in the Ysz-M 
relation. 

7. CONCLUSIONS 

We have presented the results from X-ray observations 
of a subset of 15 clusters from the first SZ selected clus- 
ter catalog from the South Pole Telescope cluster survey. 
We report the X-ray properties of this sample, includ- 
ing measurements of Mg, Tx, and Yx, and have used 
Tx and Yx to estimate the total masses of the clusters. 
We find generally good agreement between the expected 
X-ray properties of this sample and those expected from 
scaling relations. However, we find some indication of 
deviations from self-similar evolution of the Mg-T rela- 
tion compared to other local cluster samples at a level 
consistent with the explanation of an evolving gas mass 
fraction in high rcdshift clusters. 

Using the X-ray measured cluster positions and gas 
profiles, we have re-analyzed the SZ measurements to 
calculate each cluster's integrated Compton-y parameter, 
Ysz- We further use the X-ray measured gas profiles to 
deprojcct the SZ measurements so that they correspond 
to a spherical integrated Comptonization, Ysz.sph, that 
is more directly comparable to the X-ray measurements. 

We have calculated scaling relations between Ysz.sph 
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Fig. 3. — The Ysz — Yx relation for the sample. The plots show uncorrected Ygz,cyl estimates (left) and the deprojected Ygz,sph 
estimates assuming a universal temperature profile (right). The solid line shows the best fit relation. The dashed lines represent equality 
and the dotted line is the best fit relation from Arnaud ct al. (2009). SPT-CL J055f-5709 is shown as a red square and SPT-CL J0516-5430 
is shown with a blue triangle. 




Fig. 4. — The Ysz — Afsoo relation for the sample. The plots show the Ygz.sph estimates, deprojected assuming a universal temperature 
profile. Masses are estimated from the M^qq-Yx relation. The dashed line shows the best fit relation from Arnaud ct al. (2009). SPT-CL 
J0551-5709 is shown as a red square and SPT-CL J0516-5430 is shown with a blue triangle. 



and the X-ray measured quantities Yx and M^qq^Yx ■ We 
fit the Ysz,sph-Yx relation, and find a slope consistent 
with unity, 0.96 ± 0.18. Fixing this slope to 1, we re-fit 
the relation and find a normalization that implies a ra- 
tio of Ysz.sph/Yx = 0.82 ± 0.07. This normalization 
efi^ectively corresponds to the ratio between the mass 
weighted and X-ray spectroscopic weighted temperature. 
For the spherically symmetric density and temperature 
profiles assumed in this work, we would have expected a 
ratio of 0.91±0.01, consistent with predictions from other 
X-ray studies of clusters (see e.g. Vikhlinin et al. 2006; 
Arnaud et al. 2009). We therefore find a normalization 
of the Ysz,sph-Yx relation that is marginally consistent 
with, although lower than this prediction. This indicates 
that the SZ and X-ray measured pressure largely agree. 
We find, however, that the lower normalization is more 
pronounced in the z > 0.7 high-z sub-sample, with a ra- 
tio of 0.72 ± 0.14. Further X-ray and SZ observations of 
high redshift clusters are needed to determine if this is a 
real effect. 

Using the Yx measurement as a proxy for the to- 
tal cluster mass with a relation calibrated in Vikhlinin 
et al. (2009a), we find similar results when we fit a 



Ysz,sph — M^oo^Yx relation. We find a slope consistent 
with the self-similar expectation of Ysz oc M^^^ and a 
normalization consistent with the predictions from X-ray 
measurements by Arnaud et al. (2009). We have com- 
pared the Yx inferred total mass to the SZ significance- 
inferred total mass from Vanderlinde et al. (2010). Con- 
sidering only the clusters used in their cosmological anal- 
ysis, we find an average ratio of Msoo.sz to M^qq^Yx of 
0.89 ±0.06, which is within the ~ 15% systematic uncer- 
tainty on the simulation-based mass estimates in Van- 
derlinde et al. (2010). 

This work is encouraging for future studies of SZ- 
selected clusters. It suggests that SZ mass estimates are 
consistent with X-ray based mass estimates used in other 
cluster cosmology studies. This result is important for 
the use of SZ-sclcctcd cluster samples to constrain cos- 
mology, and demonstrates that the X-ray measurements 
can play a valuable role in calibrating SZ surveys. This 
work also highlights the potential power of SZ surveys to 
study cluster evolution due to the broad redshift range 
of SZ-sclected samples. The SPT has now surveyed an 
additional ^900 dcg^ to a similar depth to the ~178 dcg^ 
used in this work, and expects to cover over 2000 deg^ 
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TABLE 5 

X-RAY AND SZ SCALING RELATIONS log(y) = A+ B\og{X / Xpi-uot) + C\og(E{z)). 



Relation A B C ay ^piuot 



Lx-Tx 
Lhoi-Tx 

Lx-My,500 

Mg-Tx 

ysz,sph-yx 
ysz,aph-yx 
ysz,aph-yx 
ysz,cyi-yx 
ysz,cyi-yx 
ysz,cyi-yx 

ySZ,aph-MY,500 

Xsz.spJi-A^y.soo 



44.16 ± 0.09 


1.92 ±0.60 


1 


0.11 ±0.06 


5 keV 




44.64 ±0.14 


2.31 ±0.85 


1 


0.17 ±0.08 


5 keV 




44.63 ± 0.07 


1.16 ±0.20 


1.85 


0.11 ±0.04 


10^5 M 


13.65 ± 0.10 


1.95 ±0.66 


-1 


0.11 ±0.06 


5 keV 




14.40 ±0.08 


0.96 ±0.18 





0.09 ±0.04 


3 X 1014 


Mq keV 


14.39 ± 0.04 


1 





0.09 ±0.04 


3 X 1014 


Mq keV 


14.44 ± 0.08 


1 


-0.36 ±0.58 


0.09 ±0.04 


3 X 1014 


Mq keV 


14.50 ±0.05 


0.90 ±0.17 





0.07 ±0.05 


3 X 1014 


Mq keV 


14.49 ± 0.04 


1 





0.09 ±0.04 


3 X 1014 


Mq keV 


14.54 ±0.09 


1 


-0.35 ±0.56 


0.09 ±0.04 


3 X 1014 


Mq keV 


14.06 ± 0.10 


1.67 ±0.29 


2/3 


0.09 ±0.05 


3 X 1014 


Mq 


14.03 ± 0.04 


1.79 


2/3 


0.09 ±0.04 


3 X 1014 


Mq 



Note. — Self-similar £'(2)-scaling has been assumed here except for the Lx-^^i',500 relation, where 
the best fit evolution from (Vikhlinin et al. 2009a) is adopted, ay represents the intrinsic scatter in 
log Y. 



by December 2011. These measurements will produce 
a catalog of hundreds of massive SZ-discovered clusters 
that extend to high redshift. Joint X-ray and SZ mea- 
surements of these clusters promise to place interesting 
constraints on cluster formation and gas physics for the 
most massive young clusters in the universe. 
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APPENDIX 

NOTES ON INDIVIDUAL OBJECTS 

In this section we discuss the X-ray emission for each of the objects in the sample and note any extended substructures 
or other features in the X-ray images. Nearby radio sources, as hsted in the SUMSS survey are quoted with fluxes 
at 843 MHz, and we also discuss any particularly bright X-ray point-sources. The X-ray images (~ 15' x 15') in the 
0.5-2.0 keV band arc shown in Figures 5-19 with optical grz-images from the Magellan Baade 6.5-m telescope (see 
High et al. 2010). For SPT-CL J0546-5345 and SPT-CL J2342-5411 we instead show grz-images from the BCS survey 
since Magellan images were not available in all filters for these objects. Both X-ray and optical images are overlaid 
with SZ significance, ^, contours from VIO. The contours are from spatially filtered SPT maps where the clusters 
appear as positive significance detections, and where the spatial filtering often causes noticeable negative significance 
ringing around the cluster. The contour levels are spaced at 1.5(7 increments starting from zero with dashed contours 
for ^ < and solid contours for ^ > 0. The position of the BCG (Stalder et al., in prep.) is marked with a white cross 
in the X-ray images. Large circles and arrows show the locations of interesting X-ray features, small circles correspond 
to locations of SUMSS sources mentioned in the text. North is up and east is left. 

SPT-CL JOOOO-5748 

This cluster is the highest redshift (z = 0.74) cluster in the sample with a sharp central brightness peak (see Figure 
5). The X-ray peak is coincident with the BCG and is also associated with a SUMSS radio source with a fiux of 
40.4 ± 1.5 mJy indicative of a central AGN. The central brightness peak is broader than the Chandra PSF and we do 
not associate the central X-ray emission with AGN emission. 

SPT-CL J0509-5342 

There is a clear double-peak in the central X-ray emission of this cluster (see Figure 6). The main peak is associated 
with the central BGG. The secondary peak is located ~ 20" west of the main component (arrow) and the X-ray 
brightness ratio is 1/4. The secondary peak is also associated with two elliptical galaxies. The peak of the SZ flux is 
located between the two components. 

The main component is the most centrally compact core in this sample which indicates that a possible merger is 
in an early stage. A bright SUMSS radio source with a flux of 125.7 mJy is located 2.8' WSW of the cluster and is 
also coincident with an X-ray point source. This is not likely associated with the cluster. A possible third cluster 
component shows up as an extended X-ray source 3.5' NNW of the cluster (large circle) and is coincident with two 
elliptical galaxies in the optical image. 

SPT-CL J0516-5430 

This apparently merging cluster is very elongated along the N-S direction and also features a string of galaxies with 
the same alignment (see arrow. Figure 7). The BCG is also offset by 30" from both the X-ray and SZ peaks. There 
is a second extended X-ray component 9.5' SW of the main cluster (SW arrow). This possible subcluster also shows 
up in the SZ map as a S/N> 2 detection. The subcluster is outside of the fleld of the Magellan observation and no 
galaxy appears to be associated with this component in a Digitized Sky Survey image. 

SPT-CL J0528-5300 

This is a faint cluster (La'(0.5 — 2.0kev) — 1.9 ± 0.2 lO^'' erg s~^) without much structure and the BCG is offset by 
15" from the X-ray peak (Figure 8). The BCG is coincident with a SUMSS radio source with a ffux of 61.2 ± 2.0 mJy 
likely associated with a central AGN. 

SPT-CL J 0533- 5005 

This is another faint cluster (Lx(0.5 - 2.0kev) = 1.2 ± 0.3 10'''' erg s"^) with no weU defined core. The BCG is 
offset 45" from the central X-ray emission indicating an unrelaxed state (Figure 9). The SZ and X-ray peaks are also 
offset by 30". A local galaxy at z = 0.0147 is located 3' NW of the cluster (circle) but is not likely to affect X-ray or 
SZ measurements. 

SPT-CL JO546-5345 

The X-ray image in Figure 10 shows a substructure extending 1' SW from the main cluster (arrow) suggesting 
that a minor merger may be taking place. This elongation is aligned with an apparent extension of the SZ signal 
further supporting this scenario. One of the X-ray point sources, 6' NW, is associated with a SUMSS radio source 
with a radio flux of 19.1 ± 0.9 mJy. 

SPT-CL J0551-5709 

This apparently merging cluster has a disturbed X-ray morphology and the X-ray and SZ peaks are offset by 30". 
A SUMSS radio source with flux of 22.7 ± 1.6 mJy is located 15" east of the BCG in a region with an overdensity 
of cluster galaxies (Figure 11). Another radio source is located 5' NE of the cluster and is associated with an X-ray 
source and a local galaxy. 
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The X-ray morphology and temperature structure of this object suggest that it is projected on top of a low redshift 
structure. Extended emission can be seen both towards the S and NW, covering a large area beyond r5oo- The 
temperature decreases from 4.4 keV to 3.3 keV when increasing the aperture radius from 0.5 rsoo to rsoo, implying 
that this extended emission has very low temperature. The emission is likely associated with Abell S0552, located at 
this position. High et al. (2010) find a strong red sequence at z = 0.09 at the cluster location. The X-ray emission 
associated with this local cluster will likely bias any measurement of the temperature due to the contribution of the 
low-z flux at large radii. We attempt to correct for this bias as described in §3.5. 

SPT-CL J0559-5249 

The peak of the X-ray emission is offset to the south by ~ 25" with respect to the BCG location and SZ peak in this 
cluster (Figure 12). The X-ray emission south of the BCG is associated with a group of galaxies at this location. The 
BCG is coincident with a 50.1 ± 2.7 mJy SUMSS radio source. On larger scales, the X-ray emission and SZ signal are 
extended to the SW, clearly visible 2.5' from the cluster core (arrow in Figure 12). This cluster is clearly in a merging 
state. 

SPT-CL J2331-5051 

As noted in VIO, this cluster appears to be a pair of merging clusters with the smaller component located 2.5' SE of 
the main cluster. A spectral analysis of the smaller component indicates that it has a luminosity of 2 x lO'*'^ erg s~^, 
roughly 5% of the main cluster. As can be seen in Figure 13, the location of the X-ray emission and SZ decrement 
at the subcluster show very good agreement. This subcluster was also confirmed (see High et al. 2010) to be at the 
same redshift as the main cluster. Additionally, this cluster has a well defined X-ray cavity 18" S of the X-ray peak 
(arrow). This cluster also hosts a cool core, coincident with the BCG and a 13.3 ± 1.0 mJy SUMSS radio source. 

SPT-CL J2332-5358 

This nearby (z = 0.32) regular cluster shows no sign of significant structure in our XMM-'Newton observation with 
7 ks of good data. The location of the BCG, the X-ray peak and the peak of the SZ-signal show good agreement and 
the cluster is likely relaxed (see Figure 14). This cluster is also the only candidate in VIO with a point source detected 
in the SPT 220 GHz maps at > 5ct. We correct the SZ fiux and position of the cluster by subtracting this source from 
the 150 GHz assuming a dust spectral index as described in §5.4. 

SPT-CL J2337-5942 

This cluster shows a sharp X-ray surface brightness discontinuity NE of the cluster center (arrow in Figure 15). 
Interestingly, the BCG is also located in this region, providing tentative evidence that this is a compact gas core 
moving through the ambient cluster medium, driving a shock front. The cluster shows an overall irregular morphology 
without a clear peak in the X-ray emission. This is the cluster with the highest SZ signal in the sample and is the 
second most massive with Msoo.y^ = 7.43 ± 0.80 x 10"Mq. There are two SUMSS radio sources with 7.1 ± 1.0 mJy 
and 7.8 ±1.0 mJy fluxes 33" and 23" respectively from the cluster center. 

SPT-CL J2341-5119 

In this cluster there are no signs of significant structure within the cluster in the X-ray data (see Figure 16). 
However, 2.5' NNW of the cluster center (large circle) there is a faint extended X-ray source at the location of a bright 
galaxy indicating that this is a satellite galaxy group about to merge with the main cluster. The BCG shows good 
correspondence with the X-ray peak and the main cluster appears relaxed. There is a 7.9 ±1.1 mJy SUMSS radio 
source 6" N of the BCG which could be associated with a central AGN. 

SPT-CL J2342-54II 

This high redshift (z = 1.08) cluster shows no signs of significant merger activity and is likely relaxed. However 
there is a possible secondary component or tail to the SW, 15" from the core. 

SPT-CL J2343-5521 

In a 70 ks Chandra observation, no X-ray source was found at the location of this SPT detection. This field was also 
imaged with BCS and Magellan with no optical counterpart found. The optical data suggested that this cluster would 
have to be at 2; > 1.2 to not find an optical counterpart in either observation (High et al. 2010). The long Chandra 
observation allows us to put strict upper limits on the luminosity. Assuming the cluster is at z = 1, the 5a upper 
limit is Ljf (0.5 — 2.0kcV) < 3 x lO^^ergs s~^. The same limit assuming z = 2 is 10^"^ ergs s~^. Additionally, no cluster 
counterparts could be identified in a 3.6 /im Spitzer/IRAC observation. The X-ray, IR, and optical data all strongly 
support that this is a single band false detection by SPT. 

SPT-CL J2355-5056 

This cluster shows good agreement between the position of the BCG and the X-ray peak and shows no significant 
structure in the X-ray image (Figure 18). The SZ peak is offset ~ 35" to the south of the BCG while there is no 
evidence for extended X-ray emission in this direction. The X-ray emission is peaked in the center, likely due to a cool 
core. 
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SPT-CL J2359-5009 

The X-ray image in Figure 19 reveals a 45" offset between the location of the SPT SZ detection and the X-ray 
peak. The BCG position is offset from the SZ peak by 33". Located east of the chister is a local {z = 0.029) galaxy 
pair (circle, 7.5' E) which is clearly seen in both the X-ray and optical images. Another local galaxy [z = 0.047) is 
associated with a faint X-ray source (1.8' W), this is excluded from our X-ray analysis. A very bright point source is 
also visible in the X-ray image (3.5' NW), likely associated with an AGN. We do not find a counterpart to this source 
in the NASA/IPAC cxtragalactic database (NED) and estimate a flux of 1.6 x lO^^'^ erg s~^ cm~^ in the 0.5-2.0 keV 
band. A SUMSS radio source is located 15" SE of the BCG with a flux of 21 ± 1.1 mJy. The X-ray morphology is 
elongated along the east-west direction, possibly indicating merging activity. This is also supported by the ~ 70 kpc 
offset between the X-ray peak and the BCG position. 




Fig. 5.— SPT-CL JOOOO-5748, z=0.74 




Fig. 6.— SPT-CL J0509-5342, z=0.4626 
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Fig. 8.— SPT-CL J0528-5300, z=0.7648 




Fig. 10.— SPT-CL J0546-5345, z=1.0665 
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Fig. 12.— SPT-CL J0559-5249, z=0.6112 




Fig. 14.— SPT-CL J2332-5358, z=0.32 
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Fig. 15.— SPT-CL J2337-5942, z=0.7814 




Fig. 16.— SPT-CL J2341-5119, z=0.9983 




Fig. 18.— SPT-CL J2355-5056, z=0.35 
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Fig. 19.— SPT-CL J2359-5009, z=0.76 



